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This  dissertation  deals  with  reliable  and  accurate  extraction 

of  parameters  applicable  to  bipolar  device  and  circuit  design.  The 

extraction  is  based  on  new  methods  developed  as  a part  of  this 

study. 

Series  resistances  of  the  emitter,  base,  and  collector  are  very 
important  in  determining  the  static  and  dynamic  performance  of 
bipolar  transistors.  An  accurate  small-signal  impedance  method  for 
measurement  of  the  base  and  emitter  resistances  and  the  junction 
capacitances  is  described.  The  method  also  yields  the  current 
dependence  of  the  base  resistance.  The  method  is  demonstrated  on 
advanced  small-size  transistors  with  polysilicon  contacts.  The 
accuracy  of  the  results  is  checked  through  an  internal  self- 
consistency  test  based  on  a unique  subtraction  method  for  measuring 
the  junction  capacitances. 

A new  method  for  measurement  of  the  collector  resistance  was 
developed.  The  method  uses  the  parasitic  p-n-p  transistor  associated 
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with  the  n-p-n  integrated  circuit  transistor.  The  measured  collector 
resistance  shows  a sharp  transition  from  the  large  low-current  value 
to  the  low  high-current  value.  The  transition  is  due  to  the  high 
injection  phenomena.  The  measured  results  agree  very  well  with 
numerical  simulations.  The  method  also  yields  the  emitter  resistance 
by  combination  with  other  DC  measurement. 

The  Auger  coefficient  in  the  electron-hole  plasma  in  the 
density  range  from  1017  to  5 x 1017  cm'3  is  measured  from  the 
dependence  of  the  electroluminescence  of  the  p+-p~-n+  structure  on 
current.  The  new  method  developed  for  this  purpose  also  uniquely 
determines  the  plasma  density  as  a function  of  current.  The  range  of 
plasma  densities  studied  corresponds  to  the  electron  and  hole 
concentrations  in  the  collector  of  the  bipolar  transistor  at  high 
currents.  The  new  measured  value  for  the  Auger  coefficient  is  1.77  x 
10  3®  cm^s  which  is  about  five  times  larger  than  the  currently 
accepted  value. 
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CHAPTER  ONE 
INTRODUCTION 

As  the  complexity  of  integrated  circuits  approaches  the  level 
of  VLSI  (Very  Large  Scale  Integration),  it  is  very  important  to 
predict  the  performance  of  individual  semiconductor  devices  and 
circuits.  This  will  reduce  the  manufacturing  cost  because  correcting 
errors  in  VLSI  circuits  design  requires  a lot  of  turn-around  time. 
To  predict  the  performance  of  VLSI,  we  need  good  device  and  circuit 
simulators.  The  accuracy  of  simulators  depends  on  accurate  inputs, 
which  include  the  circuit,  physical  device  structure,  and  parameters 
of  semiconductor  devices.  Thus,  extraction  of  the  parameters  of  a 
semiconductor  device  plays  an  important  role  in  reducing  the 
development  cost  of  VLSI  circuits. 

There  are  two  kinds  of  parameters  that  describe  the  performance 
of  a semiconductor  device:  the  device  parameters  (e.g.  series 
resistances,  capacitances,  inductances)  for  the  circuit  analysis  and 
the  material  parameters  (e.g.  mobility,  bandgap,  lifetime,  etc.)  for 
the  device  analysis.  In  principle,  we  need  only  the  material 
parameters  and  the  information  about  the  structure  (geometry  and 
doping  profiles)  in  order  to  predict  the  performance  of  the  device. 
However,  there  are  still  large  uncertainties  about  the  values  for 
the  material  parameters  and  the  device  structure.  Even  if  these 
problems  were  solved,  the  simulators,  which  receive  the  material 
parameters  as  the  inputs,  cannot  give  us  an  accurate  description  of 
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a device  due  to  incomplete  understanding  of  some  physical  phenomena, 
and  also  due  to  the  multidimensional  effects  in  advanced 
submicrometer  devices. 

Considerable  effort  is  being  made  to  characterize  the  physical 
phenomena,  measure  the  material  parameters,  and  improve  the  device 
simulators,  in  order  to  extract  the  device  parameters  directly  from 
the  device  simulation  results.  However,  there  still  are 
discrepancies  between  simulations  and  experimental  results  because 
of  the  reasons  mentioned  above.  Another  approach  is  to  extract  the 
device  parameters  from  the  fabricated  devices,  which  will  serve  as 
the  input  parameters  for  the  circuit  simulator,  such  as  SPICE.  The 
second  approach  is  more  direct  and  accurate  than  the  previous 
approach,  which  starts  from  the  material  parameters  and  device 
structure  as  the  input  parameters  for  the  device  and  circuit 
simulators . 

The  direct  measurement  of  the  device  parameters  offers  a number 
of  advantages.  It  provides  information  about  devices  for  circuit 
designers  who  are  not  familiar  with  the  details  of  device  physics. 
It  may  also  help  device  physicists  to  improve  the  present 
understanding  of  the  device  operation  or,  perhaps,  to  find  unknown 
phenomena  which  are  not  included  in  the  present  models  used  by  the 
device  simulators.  The  results  of  direct  measurement  can  also  be 
used  as  criteria  to  determine  the  yield  and  uniformity  of  a certain 
device  process. 

There  are  many  techniques  for  parameter  extraction  for  bipolar 
transistor.  The  problems  involved  in  the  existing  techniques  are 
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mainly  poor  accuracy,  inapplicability  for  all  bias  levels  and 
devices,  poor  physical  meaning,  and  complexity  of  the  measurement. 
In  this  thesis,  we  suggest  new  techniques  for  measuring  the  device 
and  material  parameters  of  bipolar  transistors.  Seven  important 
device  parameters  of  a bipolar  transistor  (base  resistance,  emitter 
resistance,  collector  resistance,  collector-base  junction 
capacitance,  base-emitter  junction  capacitance,  device  temperature, 
and  junction  saturation  currents)  will  be  discussed  and  new 
measurement  techniques  for  the  above  parameters  will  be  suggested 
and  demonstrated  on  state-of-the-art  devices.  There  are  other 
parameters,  not  discussed  in  this  thesis,  that  are  required  to  fully 
describe  a transistor  model,  such  as  the  hybrid-rt  model.  The 
parameters  discussed  here,  however,  are  generally  the  most  difficult 
to  measure  and  are  the  key  parameters  that  determine  the  transistor 
characteristics  [1], 

High  speed  bipolar  transistors  operate,  for  a minimum  gate 
delay,  at  current  density  near  or  even  exceeding  the  onset  of  the 
so-called  high-current  phenomena  [2].  The  high-current  effects  in 
these  devices  occur  predominantly  in  the  lightly  doped  epitaxial 
collector  region.  If  the  transistor  operates  in  the  high-current 
regime,  an  electron-hole  plasma  can  be  generated  in  the  lightly 
doped  collector  region  [3].  The  material  parameters  in  the  electron- 
hole  plasma  may  deviate  from  their  typical  low-current  values  due 
to  carrier-carrier  scattering  and  other  effects  [4-5].  Bandgap 
narrowing,  ambipolar  mobility,  and  Auger  recombination  coefficients 
in  the  plasma  state  are  important  for  characterization  of  the  high- 
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current  regime  of  the  bipolar  transistor.  In  this  thesis,  Auger 
recombination  will  be  discussed  and  a new  measurement  technique  for 
this  parameter  will  be  introduced. 

In  Chapter  Two,  several  published  measurement  techniques  for 
the  base  resistance  will  be  briefly  reviewed.  Then,  our  new 
technique  for  measurement  of  the  base  resistance  and  the  junction 
capacitances  will  be  introduced.  Emphasis  is  placed  on  reliable 
measurements  of  small-size  transistors  with  examples  shown  for  Si 
bipolar  transistors  with  emitters  smaller  than  10  pm2.  Advantages 
over  measurements  on  large-area  test  devices  are  discussed.  The 
method  serves  as  the  basis  for  the  measurement  of  other  parameters, 
such  as  the  current-dependence  of  base  resistance,  emitter 
resistance,  space-charge  region  capacitances  of  the  base-emitter 
and  base-collector  junctions,  junction  temperature,  and  base  and 
collector  saturation  currents.  Applications  to  heterojunction 
bipolar  transistors  are  also  demonstrated. 

In  Chapter  Three,  we  will  exploit  the  contact  crowding  [6-7] 
and  emitter  crowding  [8],  to  develop  a new  model  for  an  effective 
current-dependent  emitter  resistance.  The  emitter  resistance  is 
considered  in  all  present  models  as  current  independent.  However,  we 
will  show  that  the  effective  emitter  resistance  may  vary  with 
current  and  depends  on  many  variables.  The  effective  emitter 
resistance  is  investigated  in  our  model  as  a function  of  bias, 
geometry,  sheet  resistances  of  the  base  and  emitter,  and  emitter 
specific  contact  resistance.  The  model  is  experimentally  verified  by 


measurements  on  two  extreme  device  cases. 
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In  Chapter  Four,  a new  measurement  technique  for  collector 
resistance  will  be  introduced.  The  new  technique  uses  the  parasitic 
p-n-p  transistor  of  the  main  n-p-n  transistor.  By  operating  the 
parasitic  p-n-p  transistor  in  the  upward  and  downward  modes,  we  can 
separate  various  components  of  collector  resistance.  Through  the 
accurate  value  of  collector  resistance,  we  can  also  obtain  the 
emitter  resistance.  Our  new  technique  is  demonstrated  on  various 
transistors  and  compared  with  numerical  simulations. 

In  Chapter  Five,  we  present  a new  technique  for  measuring  Auger 
coefficients  in  the  electron-hole  plasma.  Auger  coefficients  are 
determined  by  comparing  the  measured  electroluminescence  with  the 
numerically  calculated  plasma  density  of  a p+-p"-n+  diode.  The 
effects  of  series  resistance  are  removed  by  plotting 
electroluminescence  versus  current.  Another  advantage  of  the 
electroluminescence  method  is  that  it  uniquely  determines  the  plasma 
density  as  a function  of  current. 


CHAPTER  TWO 

MEASUREMENT  OF  THE  BASE  RESISTANCE 


2.1  Introduction 

The  base  resistance  rb  depends  on  the  geometry  and  the  base 
sheet  resistance  RSB , which  is  inversely  proportional  to  the 
integrated  neutral  base  impurity  density.  In  high-performance 
transistors,  the  base  width  is  very  narrow  (=  1000  A).  For  a typical 
base-emitter  junction  doping  gradient  of  1024  cnr4,  the  space-charge 
region  width  shrinks  from  about  300  & at  zero  bias  to  about  150  K at 
0.75  V [9].  The  narrowing  of  the  space-charge  region  width  widens 
the  quasi-neutral  base  and,  thus,  reduces  rb  [10].  Other  effects 
that  lower  rb  at  high  current  are  conductivity  modulation  and  base 
push-out  [H-13].  Numerical  calculation  of  the  current  dependence  of 
rb  is  possible  [14-16]  but  requires  a detailed  knowledge  of  the  base 
doping  profile,  doping  dependence  of  the  majority-carrier  mobility 
in  the  base,  and  in  small  devices,  also  a consideration  of  two-  or 
three-dimensional  effects.  The  emitter  resistance  consists  of  the 
resistance  of  the  emitter  bulk  and  the  contact  resistance.  The 
contact  resistance  is  usually  dominant  for  heavily  doped  emitter. 
The  emitter  bulk  resistance  can  be  calculated  numerically.  The 
contact  resistance  is  technology  dependent  and  has  to  be  determined 
by  measurement.  Thus,  a direct  measurement  of  the  base  and  emitter 
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resistances  is  important  not  only  for  circuit  design,  but  also  for 
verification  of  the  accuracy  of  the  process  and  device  models. 

In  this  chapter,  we  present  and  demonstrate  new  methods  for 
accurate  and  reliable  measurement  of  the  base  (rb)  and  the  emitter 
(re)  series  resistances  and  their  current  dependencies.  We  also 
determine  the  base-emitter  and  base-collector  space-charge  region 
capacitances,  temperature  of  the  device  under  test,  and  the  base  and 
collector  saturation  currents.  A self-consistency  check  (via 
measured  capacitances)  to  assure  the  accuracy  and  reliability  of  the 
measurements  on  small-size  transistors  is  described.  The  new  method 
eliminates  the  need  for  a dummy  device  normally  required  in  standard 
direct  capacitance  measurements.  An  analysis  of  the  IC-VBE  curve  to 
give  accurate  junction  temperature  and  saturation  current  is 
proposed.  The  analysis  is  based  on  a correction  of  the  measured  Ic- 
VBE  curve  for  the  effects  of  basewidth  modulation  by  VBE. 

In  Section  2.2,  we  briefly  review  conventional  methods  for 
measuring  rb  and  point  out  their  deficiencies.  In  Section  2.3,  we 
present  the  theoretical  background  for  the  new  measurement.  The 
experimental  results  for  Si  transistors  with  various  sizes  and 
application  for  GaAlAs/GaAs  heterojunction  transistor  are  also  shown 
in  Section  2.4.  Accuracy  considerations  for  measurement  of  series 
resistances  are  discussed  in  Section  2.5. 
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2.2  Review  of  Conventional  Methods 

Traditionally,  rb  has  been  a difficult  parameter  to  measure. 
The  main  reason  is  that  rb  is  modeled  as  a lumped  resistance, 
although  it  is  actually  a distributed  resistance.  Modeling  of  rb  by 
a single  lumped  resistance  simplifies  circuit  analysis  of  complex 
circuits.  As  a result  of  the  use  of  the  lumped  rb,  the  value 
obtained  for  rb  depends  strongly  on  the  measurement  technique  used 
and  on  the  transistor  operating  conditions.  Some  of  the  conventional 
techniques  are  reviewed  below. 

In  the  DC  technique  the  voltage  deviation  of  the  base  current 
from  an  exponential  ideal  dependence  is  measured.  The  voltage 
divided  by  the  base  current  gives  rb  [17].  This  technique  is  not 
valid  for  the  scaled-down  advanced  bipolar  transistors  with  poly- 
silicon contacts  because  the  voltage  deviation  includes  the  lEre 
contribution  in  addition  to  the  IBrb  term.  The  IEre  drop  may  be 
negligible  for  transistors  with  large  emitter  metal-contact  area, 
but  is  not  negligible  for  advanced  bipolar  transistors  where  re  can 
be  large  because  of  small  emitter  size.  Even  for  large  emitters  with 
metal-contacts,  there  is  a problem  with  this  technique  because  re 
increases  due  to  the  current  crowding.  In  Chapter  Three  the  increase 
in  re  will  be  discussed. 

Since  the  noise  performance  is  one  of  the  device 
characteristics  that  is  significantly  influenced  by  the  value  of  rb, 
it  is  natural  to  attempt  to  measure  noise  to  determine  rb  [18-19]. 
This  type  of  measurement  is  not  popular  because  of  problems  with 
shielding  and  noise  measurement. 
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The  semi-circle  method  uses  the  input  impedance  of  the  common- 
emitter  configuration  [20-21].  The  impedance  of  the  base-emitter 
junction  modeled  by  the  hybrid-n  transistor  model  is  plotted  on  the 
complex  impedance  plane  with  frequency  as  a variable.  The  locus 
forms  a semi-circle.  The  right  intercept  of  the  semi-circle  on  the 
real  axis  occurs  at  zero  frequency.  The  impedance  value  at  this 
point  is  the  sum  of  r^,  rb,  and  (l+3)re.  The  extrapolated  value  of 
the  left  intercept,  which  occurs  at  the  infinite  frequency,  is  the 
sum  of  rb  and  re.  The  diameter  of  the  circle  depends  on  current.  It 
is  very  large  at  small  currents  and  very  small  at  large  currents. 
The  locus,  calculated  from  a simple  one-lump  model,  is  not  an  exact 
semi-circle  [22];  thus,  a semi-circle  fit  to  the  measured  data 
yields  errors  for  intercepts.  At  high  frequencies,  the  parasitic  and 
stray  capacitances  cause  additional  distortions.  For  these  reasons, 
the  accuracy  of  this  method  is  questionable,  particularly  at  small 
currents.  Another  disadvantage  of  this  technique  is  the  time 
required  to  obtain  the  data  for  various  biases  and  the  difficulty  in 
predetermining  the  collector  current  at  which  rb  can  be  accurately 
measured. 

The  phase  cancellation  method  [23]  uses  the  input  impedance  of 
the  common-base  configuration.  The  collector  current  is  varied  until 
a certain  value  (Icm)  at  which  the  input  capacitance  is  zero.  At 
this  condition,  the  resistance  reading  is  the  sum  of  rb,  re,  and 
kT/ (qlcm) • This  method  is  not  applicable  for  transistors  with  large 
base  resistances  and  low  value  of  (3,  and  for  devices  operating  at 
high  current,  since  this  technique  assumes  rb  <<  r^.  There  are  other 
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techniques  described  in  [24],  but  these  techniques  are  not  commonly 
used  due  to  their  complexity  and  poor  accuracy.  One  approach,  which 
is  not  a direct  measurement  method,  is  numerical  simulation  [14-16, 
25].  As  mentioned  in  Chapter  One,  this  approach  requires  accurate 
knowledge  of  material  parameters  and  other  information  for  good 
accuracy.  Thus  this  approach  cannot  yield  an  accurate  absolute  value 
for  rb,  but  it  can  predict  the  general  dependence  of  the  base 
resistance  as  a function  of  current  or  voltage. 

2 . 3 Theoretical  Background 

Neugroschel  [26]  demonstrated  a small-signal  impedance  method 
to  measure  low-current  values  of  rb  and  emitter  resistance  re.  The 
method  was  demonstrated  on  relatively  large  transistors  with  the 
base-emitter  area  AgE  = 150  pm^.  in  this  chapter,  we  present  a 
modified  method  [27],  which  removes  the  uncertainties  of  the 
original  technique  and  extends  the  application  of  the  method  [26]  to 
small-size  device,  down  to  AgE  = 1 um^.  Measurement  of  rb  on  actual 
small-size  devices  is  essential  to  reveal  multi-dimensional  effects. 
The  standard  approach  of  using  large-area  devices  to  measure  rb  and 
then  do  scale-down  this  value  to  small-size  device  is  not  accurate, 
since  the  scaling  does  not  necessarily  applies  to  rb. 

The  method  presented  [26]  is  based  on  measurement  of  two  input 
impedances  of  a transistor  at  small  VBE  < 0.4  V;  the  common  emitter 
(CE)  impedance  Z^ne  and  the  common  base  (CB)  impedance  Zinb. 

To  derive  the  input  impedances  Zine  and  Zinb,  we  start  from  the 
hybrid-n  representation  of  a bipolar  transistor.  We  assume,  as 
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Neugroschel  did  [26],  that  the  effective  base-emitter  dynamic 
resistance  rn  >>  l/wCn.  This  condition  can  be  easily  satisfied  for 
any  transistor  for  VBE  < 0.4  V at  frequency  f = 50  MHz.  Other 
assumptions  used  in  the  derivation  are  given  and  justified  in 
Appendix  A.  The  expression  for  Zine  = Rlne  _ jxine  is  derived  in 
Appendix  A (see  (A. 2))  and  is  given  here  for  convenience 


Ri 


ne  = rbx  + rbi 

r cn + s j 

2 

( r 
uTt 

^ Cjne  J 

+ re 

' Cine' 

2 


(2.1) 


1/xine  ~ “Cine  - «(  Cn  + + C sc  ) = “(  Crt  + Cjc  ),  (2.2) 

where  Cjne  is  the  total  measured  capacitance  and  equal  to  the  sum  of 
all  capacitances  of  the  equivalent  circuit,  and  Cjc  is  the  total 
base-collector  space-charge  region  capacitance. 

For  the  CB  configuration  we  make  use  of  Fig.  B.l  (given  in 
Appendix  B)  and  obtain  the  following  expression  for  Z-j^g  (see  (B.2)) 

Rinb  = rbx  + rbi  + re  (2.3) 

1/Xinb  = “Cinb  = “C^.  (2.4) 


For  a typical  transistor  with  Rgg  = 10  kS/square,  rg  = rgx  + r^  >> 
re  for  small  bias  levels  and 


Rinb  = rb-  (2.5) 

There  are  two  main  differences  between  the  CE  and  CB 
expressions.  First,  Rj,ng  does  not  include  capacitances,  which  is  an 
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advantage  in  a comparison  with  Rine.  Second,  Cinb  gives  directly 
the  equivalent  base-emitter  capacitance  Cn,  while  Cine  includes  all 
capacitances  of  the  equivalent  circuit. 

2.4  Experimental  Results 
2.4.1  Base  and  Emitter  Resistances 

The  transistors  used  to  demonstrate  the  new  methods  in  this 
work  were  fabricated  using  a single  layer  poly-silicon  process  with 
a poly-silicon  emitter  contact  and  a metal  base  contact.  The 
emitter  size  in  these  devices  is  very  close  to  the  mask  size,  except 
for  a very  small  lateral  emitter  diffusion.  The  transistors  have 
non-walled  emitters  with  mask  sizes  of  2 x 4 urn2,  2 x 12  pm2,  and  2 
x 100  pm  , and  a single-stripe  base  contact.  The  devices  measured 
were  mounted  in  T0-5  headers. 

In  Fig.  2.1,  the  solid  line  shows  the  measured  dependence  of 
Rinb  ~ rb  on  VBE  for  the  transistor  with  2 x 12  pm2  emitter.  The 
data  were  obtained  at  a frequency  f = 38  MHz  with  VBC  = 0 V using  a 
Hewlett-Packard  HP4194A  impedance  analyzer.  The  measurement 
instrument  was  zeroed  with  a dummy  device  before  the  readings  of 
Zinb  were  taken.  The  measured  resistances  and  capacitances  were 
nearly  frequency  independent  in  the  range  from  20  to  40  MHz  as  shown 
in  Fig.  2.2.  Further  comments  on  importance  of  zeroing  of  parasitics 
using  a dummy  device  will  be  presented  in  Section  2.5.1. 

The  measured  Rinb  decreases  nearly  linearly  with  increasing 
VBE>  *n  an  agreement  with  numerical  predictions  [14-16].  The  sharp 
increase  in  Rinb  for  VBE  > 0.4  V is  due  to  the  violation  of  the 


Rlnb  (^) 
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vBE  (V) 


Fig.  2.1  Measured  common-base  input  resistance  as  a function  of 

base-emitter  voltage  Vgg.  The  dashed  line  is  the  linearly 
extraplated  dependence. 


Rinb  (KILO  OHM) 
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Fig.  2.2  Measured  common-base  input  resistance  for  emitter  area 

of  2 x 12  ym^  as  a function  of  frequencies  (Vg£  = 0.3  V). 
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assumption  rn  » l/<oCn  (see  Appendix  B).  Beyond  this  voltage,  (2.3) 
and  (2.4)  are  not  valid.  To  find  Rinb  for  higher  voltages,  we 
linearly  extrapolate  the  measured  Rinb  dependence,  as  shown  in  Fig. 
2.1  by  the  dashed  line.  It  was  verified  experimentally  [28],  by 
measuring  the  pinch  resistance  of  a transistor  base  with  two  base 
contacts,  that  the  extrapolation  is  valid  as  long  as  the  lateral 
voltage  drop  rbIb  « kT/q  (=  26mV  at  300  K).  The  validity  of  the 
extrapolation  was  also  verified  numerically  [14-16].  For  rbIb  > 
kT/2q,  the  current  crowding  becomes  significant,  which  causes  a 
sharp  decrease  in  rb  [29].  The  transistor  of  Fig.  2.1,  rbIb  < 5 mV 
for  VBE  < 0.9  V. 

Figure  2.3  shows  the  measured  I-V  curves  for  the  2 x 12  pm^ 
transistor.  The  base  current  IB  deviates  for  VBE  > 0.8  V from  its 
nearly  ideal  exp(qVBE/kT)  dependence  by 

AV  = rbIB  + reIE.  (2.6) 

Equation  (2.6)  is  valid  for  any  VBE,  provided  that  IB  is  due  to  the 
recombination  in  the  heavily  doped  emitter,  which  remains  in  low 
injection.  The  effect  of  re  on  AV  is  not  negligible,  even  if  re  « 
rbi,  since  IE  = (l+e)reIB. 

To  find  re,  we  combine  AV  from  Fig.  2.3  with  rh  = Rinb  from 
Fig.  2.1.  For  the  VBE  range  from  0.8  V to  0.85  V,  we  find  that  AV 
varies  from  7 to  25  mV.  Using  (2.6)  we  find  that  re  = 10  S.  A more 
accurate  SPICE  analysis  (see  Chapter  Three),  using  a filament 
transistor  model  with  20  sections,  yields  only  a small  correction 
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vBE(V) 


Fig.  2.3  Measured  dependencies  of  the  base  and  collector  currents  on 
VfiE  f°r  two  values  of  Vg£.  The  dashed  line  extending  the  Ig 
plot  is  the  ideal  base  current  proportional  to  exp(qV/kT). 
The  dotted  curves  are  for  the  substrate  current  Is  for  Vcs 
= 0 V (upper  curve)  and  0.2  V (lower  curve). 
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shown  in  Fig.  2.1  by  the  dashed-dotted  line  and  gives  re  = 10.2  J3. 
Similar  analysis  for  a 2 x 100  pm2  transistor  gave  re  = 2.2  52 
Assuming  that  re  is  due  mainly  to  the  contact,  this  corresponds  to 
the  specific  emitter  contact  resistance  of  pc  = reABE  = 236  S3. pm2 
and  pc  = 360  S3. pm2,  respectively. 

In  spite  of  the  fact  that  current  crowding  occurs  for  VBE  > 
0.85  V,  the  emitter  current  through  the  emitter  contact  will  still 
be  essentially  uniform  (see  Chapter  Three),  and  re  will  remain  at 
the  low-current  value  determined  above  (assuming  that  pc  is 
constant).  This  allows  calculating  rb  with  (2.6)  as  a function  of 
VBE  or  the  collector  current  Ic.  The  results  are  shown  in  Fig.  2.4 
(a)  and  (b). 

The  kink  in  the  IB  versus  VBE  curve  in  Fig.  2.3  at  VBE  = 1.1  V 
is  due  to  the  forward-biasing  of  the  base-collector  junction  by  the 
voltage  drop  ICRC  on  the  collector  resistance  Rc,  which  contributes 
to  IB  in  addition  to  the  contribution  from  the  base-emitter  junction 
(see  Chapter  Four).  Corresponding  kinks  seen  in  the  rb  versus  VBE 
and  rb  versus  Ic  plots  in  Fig.  2.4  (a)  and  (b)  are  due  to  the 
modulation  of  rb  by  injection  from  the  forward  biased  base-collector 
junction.  The  IcRc  drop  will  also  forward  bias  the  collector- 
substrate  junction  causing  Ig  to  flow  into  the  substrate.  The  kink 
in  IB-VBE  plot  can  be  also  caused  by  the  emitter  current  of  a 
parasitic  p-n-p  transistor  flowing  into  the  base.  The  substrate 
current,  Ig,  measured  for  VBC  = 0 V and  two  values  of  the  collector- 
substrate  reverse  bias,  Vcs  = 0 and  0.2  V,  is  shown  in  Fig.  2.3  by 
the  dotted  lines.  It  is  clear  from  Fig.  2.3  that  Ig  « Ig  (note 
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(a) 


Fig.  2.4  The  total  base  resistance  as  a function  of  (a)  Vru  and 

(b)  IC. 
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(b) 


Fig.  2.4  (continued) 
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also  that  IB  is  the  same  for  both  values  of  Vcs).  This  clearly 
shows  that  the  kink  in  the  Ig-Vgg  curve  is  not  caused  by  Ig  for  the 
transistor  of  Fig.  2.3.  The  fact  that  Ig  flows  into  the  substrate 
and  not  out  of  the  substrate  was  verified  by  observing  the  polarity 
of  Ig.  The  deviation  of  Iq  from  the  ideal  exponential  dependence  is 
determined  not  only  by  ohmic  losses  (AV),  but  also  by  other  high- 
current  effects.  As  a result,  no  visible  kink  in  ln(Ic)  versus  VBE 
is  observed,  although  a small  kink  is  observed  when  the  plot  of  Ic 
versus  VBE  is  displayed  on  a linear  scale. 

The  broken  lines  in  Fig.  2.4  (a)  and  (b)  show  the  rb  plots  for 
^BC  = 1 V.  Compared  to  the  zero  bias  case,  the  reverse  bias  narrows 
the  base  thickness,  which  increases  the  low-current  rb,  but 
decreases  the  current-crowding  onset  voltage.  The  low-current 
portion  of  the  rb  dependence  for  VBC  = 1 V could  not  be  measured 
with  the  present  experimental  setup.  At  high  currents,  rb  will 
probably  saturate  at  a value  close  to  rbx.  For  VBE  > 1.1  V, 
however,  our  procedure  for  measuring  rb  is  not  very  reliable.  One 
reason  is  a possible  heating  of  the  junction  due  to  large  Iq.  In 
addition  (2.6)  uses  the  assumption  that  IB  is  due  to  the 
recombination  in  the  emitter  only.  This  assumption  may  be  violated 
because  of  the  large  density  of  holes  and  electrons  recombining 
within  the  base  and  collector.  The  recombination  at  high  currents 
has  to  account  for  the  bandgap  narrowing  induced  by  the  electron- 
hole  plasma  [30],  which  further  increases  the  current. 
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Results  for  Rinb  for  the  three  transistors  measured  are  shown 
in  Table  2.1.  Since  the  emitters  in  these  transistors  are  non- 
walled , the  calculation  of  the  base  sheet  resistance  from  Rinb  ls 

not  straightforward  and  depends  on  the  layout  as  well  as  other 
variables  [25,31]. 

2-4.2.  Extrinsic  Base  Resistance  rbvr 

Because  rbx  can  be  considered  as  current  independent,  except 
for  very  high  currents,  it  can  be  conveniently  determined  from  the 
zero-bias  values  of  Rine(0)  and  Rinb(0).  If  Cy  « Cn,  which  is  a 
good  approximation  because  the  base-emitter  space-charge  region 
capacitance  per  unit  area  is  about  10  times  larger  than  that  of  the 
base-collector  junction  (see  Tables  2.2  and  2.3);  then, 

rbx  ~ Rinb(°>  - rbi(0) 

Rine(0)  ~ Rinb(0)  I + Cy)/Cine]^ 

= i ~ • (2.7) 

1 [(Cji  + Cjj)/Cine]2 

Because  of  coupling  between  the  active  and  inactive  base,  a 
separation  of  the  total  base  resistance  rb  into  rbi  and  rbx  is 
possible  only  for  transistors  with  walled  emitters  or  for 
transistors  with  very  long  emitter  stripe.  Since  transistors  with 
walled  emitters  were  not  available,  the  measurement  was  done  for  a 2 
100  pm2  transistor  and  yielded  rbx  =6  5?. 
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Table  2.1  Summary  of  the  measured  results  for  the  base  resistance 


Emitter  Mask 
Size 
(pm2 

cs 

(fF) 

RS 

(G) 

cinb 

(fF) 

Acpar 

(fF) 

^inb(O) 

(52) 

2x4 

30.4 

1023 

21 

+9.4 

2143 

2 x 12 

70.8 

904 

73.5 

-2.7 

839 

2 x 100 

568 

129.9 

589 

-21 

121 

Table  2.2  Summary  of  the  measured  base-collector  space-charge  region 
capacitances.  The  zero-bias  capacitance  Cic(0)  is  obtained 
from  (2.8).  J 


Emitter  Mask 
Size 
(pm) 

CCE 

(fF) 

CCB 

(fF) 

Cjc(0) 

(fF) 

aBC 

(pm^) 

2x4 

117 

95 

23 

8 x 10 

2 x 12 

201.5 

143.5 

58 

16  x 10 

2 x 100 

1048 

694 

354 

104  x 10 

Table  2.3  Summary  of  the  measured  base-emitter  space-charge 
region  capacitances.  and  C^g  are  the  measured 
capacitances  for  the  inverted  emitter-down  transistors. 


Emitter  Mask 
Size 

CCE 

CCB 

Cjg(O) 

(Um2) 

(fF) 

(fF) 

(fF) 

2x4 

117 

96 

21 

2 x 12 

201.5 

128 

73.5 

2 x 100 

1059 

470 

589 
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2'4'3*  Base-Emitter  and  Base-Collector  Space-Charge  Capacitance 

To  determine  the  base-collector  space-charge  region  (junction) 
capacitances  Cjc,  ve  measure  the  input  CE  and  CB  capacitances.  The 
measured  input  capacitances  will,  in  addition  to  the  intrinsic 
capacitances  Cine  and  Cinb  given  by  (2.2)  and  (2.4),  include  the 
parasitic  capacitance.  The  parasitic  capacitance  Cp,  due  to  the 
bondpads  and  metallization  over  the  oxide,  will  appear  in  the 
equivalent  diagrams  in  Fig.  2.5  between  the  terminals  and  the 
substrate  which  is  left  floating.  It  can  be  shown  that,  for 
conditions  in  which  (2.3)  and  (2.4)  are  valid,  the  parasitic 
capacitances  Ccs,  including  the  collector-substrate  capacitance,  can 
be  coupled  into  one  effective  capacitance  Cpar  between  the  input 
base  and  emitter  terminals  for  both  the  CE  and  CB  configurations  as 
shown  in  Fig.  2.5  (a)  and  (b).  The  total  input  capacitances  are  then 

CCE  = Cine  + cpar  and  CCB  = cinb  + Cpar.  Thus,  using  (2.2)  and 
(2.4)  we  obtain 

°CE  " CCB  = (Cine  + cpar)  ~ (Cinb  + Cpar)  = Cjc.  (2.8) 

The  base-collector  capacitance  is  obtained  simply  by  subtracting  the 
two  measured  input  capacitances.  Similarly,  the  base-emitter 
junction  capacitance  Cje  is  obtained  from  measurements  of  the 
inverted  emitter-down  transistor. 

One  important  difference  between  the  new  method  and  the 
standard  direct  measurement  of  junction  capacitances  is  that 
accurate  direct  measurement  requires  perfect  zeroing  by  a dummy 
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Fig.  2 


(a) 


(b) 

5 The  equivalent  circuits  of  common-emitter  (a)  and  common 
-base  (b)  for  a transistor.  C and  Ccs  at  the  inputs 
are  the  stray  capacitances  between  the  bonding  pads  and 
the  substrate  and  the  collector-substrate  capacitance, 
respectively. 
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device.  This  is  not  possible  for  small  devices  in  which  the 
parasitic  capacitances  are  much  larger  than  the  intrinsic 
capacitances  because  of  variations  over  the  wafer  in  oxide  thickness 
and  geometry,  and  also  because  the  metal  coverage  for  the  dummy  is 
different  from  that  for  the  active  transistor  in  the  vicinity  of 
the  contact  holes.  Additional  error  may  be  due  to  the  stray 
capacitances  of  the  probe  tips  for  on-wafer  probing.  Because  our 
method  does  not  require  the  dummy  device,  the  measurement  error  due 
to  zeroing  is  practically  eliminated.  The  measurement  can  be  done 
without  disconnecting  the  probes  by  using  a jumper  to  short  the 
collector  to  the  emitter  or  to  the  base.  Both  input  capacitances 
can  be  measured  with  an  accuracy  of  about  ±0.3  fF  allowing 
measurement  on  devices  with  ABE  = 1 pm2.  The  measuring  frequency 
can  be  as  low  as  1 MHz,  or  less,  making  on-wafer  measurements 
relatively  easy. 

The  results  for  the  zero-bias  capacitances  Cjc(0)  and  Cje(0)  in 
three  different  transistors  measured  at  a frequency  of  1 MHz  are 
summarized  in  Tables  2.2  and  2.3.  The  results  of  on-wafer 
measurements  for  MOSAIC  III  devices  are  shown  in  Table  2.4.  Column  3 
in  Table  2.4  shows  the  results  of  the  direct  measurement  method, 
Cje(Mot),  for  comparison.  The  agreement  is  excellent  except  for  the 
smallest  devices.  This  shows,  that  both  techniques  are  accurate.  The 
substraction  method  using  (2.8)  will  be,  however,  more  reliable  for 
smaller  devices  because  of  the  reasons  discussed  briefly  above.  The 
method  can  be  used  for  any  reverse  bias. 


The  applicability  for 
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Table  2.4  Summary  of  on-vafer  measurements  for  MOSAIC  III  devices. 


Emitter 
Mask  Size 
(ym2) 

Cjc 

(fF) 

Cjc (Mot) 
(fF) 

Cje 

(fF) 

Cje(Mot) 

(fF) 

1.5  x 4 

15.8 

17.0 

21.1 

19.0 

1.5  x 4 

16.2 

17.0 

20.6 

20.0 

2 x (1.5  x 4) 

30.7 

32.0 

39.4 

39.0 

1.5  x 8 

25.4 

26.0 

38.1 

38.0 

1.5  x 12 

34.0 

34.0 

57.2 

56.0 

2 x (5.75  x 16.5) 

188.8 

186.0 

392.6 

393.0 

2 x (5.75  x 16.5) 

189.3 

186.0 

396.7 

391.0 

6 x (2  x 19) 

321.7 

322.0 

633.7 

613.5 

2 x 12 

36.5 

37.0 

67.2 

66.0 

2x8 

27.9 

27.0 

44.7 

45.0 

2 x (2  x 4) 

33.1 

34.0 

47.9 

47.5 

2x4 

18.1 

18.0 

25.5 

24.0 

2x4 

17.5 

17.0 

25.9 

24.0 
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forward  bias  is  limited  by  the  requirement  wCrt  » r^,  and  breaks 
down  for  VBE  > 0.4  V (see  Fig.  2.1). 

Cje  can  be  accurately  extrapolated  to  higher  voltages  either 
numerically  [32]  or  analytically  [9]  from  the  measured  Cje(0).  The 
calculations  for  VBE  > 0.8  V should  include  the  effects  of  the 
electron-hole  plasma  in  the  SCR  [30].  This  effect  is  generally 
ignored,  but  recent  calculations  show  that  the  electron-hole  plasma 
will  increase  the  capacitance  by  as  much  as  a factor  of  three  to 
four  in  a comparison  with  standard  values  [33]. 

2‘4,4,  Measurement  of  Junction  Temperature  and  Saturation  Currents 
The  ideal  collector  current  of  a bipolar  transistor  and  the 
base  sheet  resistance  RSB  « rbi  can  be  expressed  by  standard 
relations: 

Aqnj2  exp  (qVBE/kT) 

Ic  — 

Vb  Nb(x) 

dx 

JO  Dn(x) 


Aqni2Dn  exp(qVBE/kT) 

~ = IC0exP(qvBE/kT> 

NB(x)dx 

Jo 


(2.9) 


rsb 


pb 

qup(x)NB(x) 


-1 


P>b 

NB(x)  dx 

o J 

(2.10) 
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where  NB(x)  is  the  position  dependent  base  doping,  Dn  is  the 
spatially  averaged  minori tv— carrier  diffusion  coefficient  in  the 
base,  and  Pp  is  the  spatially  averaged  majority-carrier  mobility 
in  the  base.  The  junction  temperature  can  be  obtained  using  (2.9) 
from  kT/q  = dVBE/d(lnIc) . This  is  strictly  valid  only  if  the 
width  of  the  quasi-neutral  base  Wb  remains  constant  with 
increasing  VBE.  But,  as  Fig.  2.1  clearly  indicates,  Wb  increases 
(measured  Rinb  decreases)  with  increasing  VBE . Since  Ic  depends 
inversely  on  the  integrated  neutral  base  impurity  density,  the 
measured  Ic  varies  with  VBE  ns  ICOexp(qVBE/mkT)  with  the  voltage 
dependent  slope  factor  m(VBE)  > 1.  Thus, 

dVBE/d(lnIc)  = m(VBE)kT/'i . (2.11) 

Equation  (2.11)  cannot  be  used  to  obtain  the  accurate  junction 
temperature,  because  the  voltage  dependent  slope  factor  m(VBE)  is 
not  known.  To  obtain  the  correct  temperature  from  Iq,  we  have  to 
make  a correction  to  account  foi  the  Wb  variation  with  VBE.  This 
can  be  done  approximately  by  taking  the  ratio  of  (2.9)  and  (2.10) 

Ic/rsb  = Aq2ni2  ( Dnpp  ) exp(qVBE/kT) . (2.12) 

The  parameters  Dn  and  Mp  are  somewhat  voltage  dependent,  but  it 
can  be  shown  numerically  that  in  the  VBE  range  from,  e.g.,  0.4  V to 
0.8  V,  the  dominant  factor  causing  the  Rinb  decay  shown  in  Fig. 

2.1  is  the  integrated  base  doping.  Taking  Dnpp  as  VBE  independent 
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gives  I(VRsb  - Cexp(qVBE/kT),  where  C is  a temperature- 
independent  constant,  and 

dVBE/d(ln(Ic/RSB))  = ItT/q.  (2.13) 

The  slope  given  by  (2.13)  is  independent  of  VBE  and  yields  the 
correct  junction  temperature  T.  A correction  for  the  voltage 

dependence  of  Dnpp  can  be  made  numerically,  if  the  base  doping 
is  known. 

To  illustrate  the  procedure  we  give  here  the  results  of 
analysis  for  a 2 x 100  pm2  transistor.  The  saturation  current  Ic0 
obtained  from  the  measured  IC-VBE  curve  depends  on  the  range  of 
voltage  chosen.  The  least  square  fit  for  the  range  from  0.34  V to 
0.7  V yields  kT/q  = 25.84  mV  (T  = 299.6  K)  and  Ic0  = 2.6  x l0-16  A. 
Using  the  range  from  0.5  V to  0.7  V gives  kT/q  mV  (T  = 300.76  K ) 
and  IC0  = 2.83  x 10-16  A.  The  I(;/rbl  dependence  gives  a straight 
line  with  the  slope  kT/q  = 25.51  mV  (T  = 295.94  K)  and  ICQ  = 2.2  x 
10-16  a using  rb.(0)  , rb(0)  - rbx  . m a _ 7 a . 114  a.  Using  the 
correct  T we  find  the  slope  factor  m = 1.012  - 1.016  for  the 
measured  Ic.  The  ambient  temperature  monitored  by  a precision 

thermistor  with  an  accuracy  of  ±0.2  K was  T = 296.26  K. 

The  discrepancy  of  about  4 K between  the  actual  T and  that 
obtained  from  the  Ic  dependence  corresponds  to  the  ratio  of  the 
intrinsic  carrier  concentrations  of  about  n^OOO  K)/ni2(296  K)  = 2. 
This  discrepancy  is  important  if  we  want  to  verify  the  accuracy  of 
device  simulators  or  analytical  models  by  comparing  the  measured  IfJ 
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and  IB  with  calculations  (which  require  the  temperature  as  an 
input).  The  discrepancies  between  Ico  from  the  IC-VBE  and  (Ic/rbi)- 
VBE  curves  will  be  much  larger  at  smaller  temperatures.  For 
example,  at  T = 100  K,  an  extrapolation  of  Ic  from  VBE  = 1 V to  VBE 

= 0 V with  the  slope  factor  of  m = 1.015  will  overestimate  ICQ  by  a 
factor  of  about  5. 

Saturation  currents  depend  on  material  parameters,  such  as 
mobility,  bandgap  narrowing,  doping  profile,  etc.,  and  are, 

therefore,  the  main  measured  parameters  from  which  the  above 
transport  parameters  are  extracted.  This  cannot  be  done  without 
accurate  and  reliable  determination  of  the  saturation  currents  and 
the  corresponding  device  temperature. 

The  measured  In  IB-VBE  curves  for  the  devices  tested  also  shows 
m > 1.  Thus,  a direct  extrapolation  gives  an  error  in  IBQ.  The  error 
is  minimized  by  taking  IB0  = 1^/^,  where  ICQ  is  found  from 

(2.16)  and  $„ax  is  the  maximum  DC  current  gain.  It  is  reasonable  to 

assume  that,  for  Vdi?  at  which  — ft  r.  u 

EE  ai  wnicn  p _ the  space-charge  region 

recombination  is  minimized. 

2‘2 * 4'5  gg£istance  Measurements  for  Hetrojunction  Bipolar  Transire 
Resistance  analysis  of  HBTs  using  the  DC  I-V  curves  is 

particularly  difficult  because  typical  I-V  curves  are  very  non- 

ideal. This  makes  even  an  estimate  of  resistance  using  DC  methods 

difficult.  To  illustrate  this,  we  show  in  Fig.  2.6  a typical  IB-VBE 
curve  for  an  n-p-n  GaAlAs/GaAs  HBT  that  follows  an  IB  * exp(qV/mkT) 
dependence  with  m = 1.5  to  1.9.  This  indicates  that  IB  is  dominated 
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VBE(V) 


Fig.  2.6  Iq-Vbe  and  IB-VBE  characteristics  for  a heterojunction 

GaAlAs/GaAs  transistor.  The  dashed  line  is  a reference  line 
with  an  ideal  exponential  voltage  dependence. 
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by  the  SCR  recombination  current.  The  IC-VBE  curve  gives  m = 1.05. 
The  transistor  was  fabricated  by  LPE  growth  method  and  has  a 
circular  geometry  with  a diameter  of  60  urn.  The  base  width  is  about 
0.2  ym  and  is  uniformly  doped  with  about  5 x 1018  acceptors/cm8. 
Analysis  of  series  resistances  for  transistor  with  such  non-ideal 
base  current  using  DC  techniques  is  not  possible,  since  they  require 
an  ideal  exp(qV/kT)  component  as  a reference  current.  Our  AC 

technique  does  not  have  this  restriction  and  is  ideally  suited  for 
the  analysis  of  HBTs.  The  results  are  shown  in  Fig.  2.7  (a)  and 
(b).  Figure  2.7  (a)  shows  the  measured  resistances  Rinb  and  Rine  as 
a function  of  VBE.  Figure  2.7  (b)  displays  the  resistances  rbi  and 
rbx  obtained  from  (2.1),  (2.3),  and  (2.7).  In  contrast  to  the  Si 
case  of  Fig.  2.1,  the  rbi  for  the  HBT  is  almost  current  independent 
because  the  base  is  heavily  and  uniformly  doped,  which  makes  the  RgB 
modulation  due  to  VBE  negligible.  Because  of  large  bandgaps  of 
GaAlAs  and  GaAs,  the  applicability  of  the  ac  method  extends  about 
VBE  = 1 V.  The  large  base  doping  makes  the  intrinsic  and  extrinsic 
base  resistances  comparable,  and  they  can  be  both  determined 
accurately.  The  extrinsic  resistance  rbx  includes  the  base  contact 
resistance,  which  appears  to  be  very  important  in  limiting  the  speed 
of  HBT  [34-35].  The  emitter  resistance  was  negligible  for  the  very 
large  test  transistor  and  could  not  be  found  accurately.  However, 
for  small-size  HBTs,  which  were  not  available  for  testing,  all  three 
resistances  can  be  measured. 
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Fig.  2.7  Measured  results  for  the  heterojunction  transistor. 

(a)  Measured  input  resistances  for  the  common-base  and 
common-emitter  for  the  HBT  from  Fig.  2.6. 

(b)  The  base  intrinsic  and  extrinsic  resistances  versus  Vnc- 

from  (a).  BE 
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2.5  Discussion 

2,5,1  Agcuracy  Considerations  for  Measurement  of  Base  Resistance 

The  procedure  for  measuring  the  base  resistance  using  the 
common-base  configuration  requires  compensation  of  parasitics  by 
zeroing  with  a dummy  device.  For  a perfect  zeroing,  the  measured 
values  Rs  and  Cs,  are  equal  to  the  intrinsic  parameters  Rlnb  and 
cinb  given  by  (2.3)  and  (2.4).  Since  our  method  uses  frequency  f < 
50  MHz,  the  transmission  line  effects  of  the  parasitics  can  be 
neglected  and  they  can  be  considered  to  be  nearly  pure  capacitors. 
Since  a perfect  zeroing  for  Cpar  is  not  possible,  an  error  AC_ 

p3.1T 

gives  Cs  = Cinb  + ACpar,  and 


Cinb 

2 

f Cinb  1 

■ ^inb+A  Cpar  > 

~ Rinb 

l cs  J 

To  be  sure  that  no  error  in  Rinb  is  caused  by  assuming  that  ACnar.  = 

psr 

0,  we  need  a self-consistency  test  for  the  zeroing  procedure  that 
would  give  us  the  residual  uncompensated  or  overcompensated  AC 

psr 

(either  positive  or  negative)  and  Cinb.  This  is  particularly 

important  for  small  devices  where  Cpar  » cinb.  The  correct  Rinb 
can  be  then  obtained  from  (2.14),  in  spite  of  imperfect  zeroing.  The 
compensation  offset  ACpar  is  obtained  from  two  measured  values  as 

^par  = cs  - cinb*  As  pointed  out  above,  Cinb  for  the  inverted 
emitter-down  transistor  is  found  very  accurately  from  Cjnb  = CCE  - 

CCB.  This  measurement  is  immune  to  parasitics  because  of  the 
subtraction  involved. 
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Table  2.1  shows  the  measured  values  Rs,  Cs,  and  Cinb  = Cje  (see 

Table  2.3).  Notice  that  ACpar  is  positive  for  the  2 x 4 pm2  device, 

but  is  negative  for  the  other  two  devices.  The  last  column  in  Table 

2.1  shows  the  corrected  values  for  rb  = Rinb  obtained  using  (2.14). 

We  have  found  that  most  reliable  and  reproducible  results  are 

obtained  by  measuring  both  the  dummy  device  and  the  active  device  in 

a parallel  equivalent  mode.  The  capacitance  and  conductance  of  the 

dummy  device  are  then  subtracted  from  those  for  the  active  device 

and  converted  to  the  series  values  R_  and  Co. 

s 

2,5,2  events  on  the  Open-Collector  Method  for  Measurement  of  the 
Emitter  Resistances  rc  -■  K 

The  open-collector  method  is  the  most  commonly  used  techniques 
to  measure  re.  The  method  requires  only  DC  measurement  of  the 

collector-emitter  saturation  voltage  VCES  as  a function  of  IB  with 
the  collector  open-circuited.  The  emitter  resistance  is  then 
calculated  from  re  = dVCES/dIB  [24].  Some  modifications  of  this 
procedure  are  suggested  in  [36-37].  Figure  2.8  shows  the  VCES 
dependence  on  Ifi  for  the  2 x 12  pm2  transistor.  The  reciprocal  slope 
of  the  dashed  line  taken  at  VCES  = 30  mV  gives  dVCES/dIB  = 17.5  q. 
The  dependence  is  not  linear  and  suggests  that  re  decreases  with 
increasing  IB.  For  example,  dVCES/dIB  = 14  8 at  VCES  = 80  mV.  It 
is  easy  to  show  that  neither  of  these  values  is  possible.  To 

demonstrate  that  we  go  to  Fig.  2.3  and  translate  the  measured  IB-VBE 
curve  to  the  left  by  the  voltage  drop  (-AV)  = IEre  for  VBE  < 1 v. 
The  corrected  IB,  using  the  open-collector  values  for  re,  will  be  on 
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Fig. 


VCES  (mV) 


2.8  Collector-emitter  saturation  voltage  as  a function  of  base 
current  for  the  2 x 12  ym^  for  measurement  of  emitter 
resistance  by  flyback  method. 
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the  left  of  the  ideal  dashed  line  (even  without  considering  the  rbIb 
contribution  to  AV),  which  is  not  possible.  The  open-collector 
method  clearly  overestimates  re.  The  correct  value  for  re  = 10.2  2 
was  found  previously  using  our  method.  Another  technique  for 
measuring  re  will  be  suggested  in  Chapter  four. 


2 • 6 Conclusions 

New  measurement  methods  for  accurate  analysis  of  submicrometer- 
emitter  bipolar  transistors  were  presented.  Measurement  of  the  base 
resistance  for  a wide  range  of  collector  current  was  demonstrated. 
The  new  method  for  the  measurement  of  the  space-charge  region 
capacitances  does  not  require  the  use  of  a dummy  device  for  the 
compensation  of  the  parasitics,  which  avoids  zeroing  errors 
occurring  in  standard  capacitance  measurements.  The  accuracy  of 
about  ±0.3  fF  for  the  measurement  of  the  base-emitter  and  base- 
collector  capacitances  can  be  employed  to  study  the  effects  of 
process  variations  in  order  to  optimize  the  speed.  The  ability  to 
measure  small-size  transistors,  instead  of  using  large  test  devices, 
can  be  used  to  investigate  yield,  uniformity  across  the  wafer,  and 
multi -dimensional  effects.  The  measurements  suggested  in  this  paper 
can  be  performed  at  frequencies  from  about  1 to  about  40  MHz.  This 
relatively  low  frequency  range  makes  on-wafer  testing  feasible. 


CHAPTER  THREE 

MODELING  OF  THE  EFFECTIVE  EMITTER  RESISTANCE 
3.1  Introduction 


In  the  most  widely  used  common-emitter  (CE)  configuration  the 
emitter  resistance  is  amplified  by  the  CE  current  gain  ((3).  Thus,  an 
accurate  analysis  of  the  emitter  resistance  is  necessary  for  exact 
transistor  modeling.  The  emitter  resistance  is  usually  considered  to 
be  dependent  only  on  the  physical  parameters  of  a transistor. 
However,  as  will  be  shown  later,  the  emitter  resistance  also  varies 
with  current  because  of  contact  and  emitter  crowding.  The  emitter 
resistance  at  low  current  levels  is  simply  determined  by  the  emitter 
contact  resistance.  The  emitter  resistance  at  the  operating  point 
where  the  emitter  crowding  is  significant  has  two  components.  The 
first  component  is  from  the  interface  layer  of  the  contact,  the 
second  one  is  from  the  semiconductor  layer  under  the  contact  hole. 

In  Section  3.2  we  introduce  physics  necessary  to  characterize 
the  contact  resistance.  In  Section  3.3  we  develop  new  model  for  the 
accurate  evaluation  of  the  emitter  resistance  as  a function  of 
current.  The  extraction  method  for  the  emitter  resistance  based  on 
the  new  model  is  suggested  in  Section  3.4.  Section  3.5  presents  the 
computer  simulation  results  and  the  experimental  results. 
Conclusions  will  be  given  in  Section  3.5. 
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3.2  Basic  Physics  for  the  Contact  Resistance 
3.2.1  Classification  of  Contacts 

The  contacts  of  semiconductor  devices  can  be  classified  into 
two  types,  the  horizontal-type  or  the  vertical-type,  depending  on 
the  direction  of  the  current  flow  in  the  semiconductor  under  the 
planar  contact.  For  a better  illustration,  a cross-section  of  a 
bipolar  transistor  is  shown  in  Fig.  3.1.  Near  the  base  contact 
region,  the  current  flows  mostly  horizontally,  which  causes  current 
crowding  at  the  contact  edge.  This  characterizes  the  horizontal-type 
contact.  Since  most  of  the  collector  current  in  the  integrated 
circuit  transistor  flows  through  the  heavily-doped  buried  layer,  the 
current  under  the  collector  contact  flows  almost  vertically.  Thus, 
the  collector  contact  is  an  example  of  a vertical-type  contact.  The 
emitter  contact  belongs  to  both  types,  depending  on  the  bias  level 
and  other  factors.  At  the  bias  level  where  the  emitter  crowding  is 
not  significant,  the  emitter  contact  acts  as  a vertical  contact.  As 
the  emitter  crowding  becomes  significant,  the  emitter  contact  may 
become  the  horizontal-type  contact,  or  in  between  the  horizontal  and 
vertical  types.  Most  of  the  contacts  in  MOS  belong  to  the 
horizontal-type.  Many  efforts  were  made  to  characterize  the 

horizontal-type  contact  for  MOS  and  other  devices  [38-45].  The 
characterization  for  the  emitter  contact  is  more  complicated  than 
for  a pure  horizontal-type  because  of  the 


reasons  mentioned  above. 


COLLECTOR 

CONTACT 


EMITTER 

CONTACT 


BASE 

CONTACT 


Fig.  3.1  Cross-section  of  a planar  n-p-n  transistor  with  the  main 
current  flows  indicated  by  arrows. 
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3*2.2  Specific  Contact  Resistance 

The  specific  contact  resistance  (pc)  in  52-cm2  is  defined  as  the 
resistance  of  unit  area  of  the  interfacial  layer  between  two 
conduction  layers  and  is  expressed  as  the  ratio  of  voltage  drop 
across  the  interfacial  layer  to  the  current  density.  The  carrier 
transport  through  the  interfacial  layer  is  governed  by  the 
thermionic  emission  over  the  barrier  and  tunneling  through  the 
barrier  [46].  The  relative  magnitude  of  the  two  components  depends 
on  the  barrier  height,  temperature,  doping  concentration,  and  other 
factors.  For  low  doped  samples,  the  thermionic  emission  component 
dominates  over  the  tunneling  component.  When  the  doping 
concentration  is  high,  the  depletion  region  width  near  the  interface 
in  the  semiconductor  becomes  very  narrow.  Therefore,  carriers  can 
easily  tunnel  through  the  barrier. 


3'2*3  Contact  Crowding  and  Transmission  Line  Model 

Using  the  definition  of  pc,  we  can  easily  estimate  the 
resistance  value  for  the  vertical  contact.  But,  the  contact 
resistance  for  the  horizontal  contact  cannot  be  easily  expressed  by 
Rc  = pc/Area,  because  the  current  density  through  the  contact  hole 
is  not  uniform  due  to  the  so-called  contact  crowding  [6-7].  Figure 
3.2  shows  the  cross-section  of  a contact  and  its  distributed 
equivalent  circuit.  In  Fig.  3.2  (b), 


dR  = 


V 


dx 


(3.1) 
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(a) 


Fig.  3.2  Schematic  illustration  of  a contact  hole  for  the  modeling 
of  contact  crowding. 

(a)  Cross-section  of  a contact  hole.  Current  flows  into  the 
left  edge  of  the  contact  hole  and  flows  out  from  the 
right  edge. 

(b)  Transmission-Line  Model  for  the  contact  depicted  in 
(a). 
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(b) 


Fig.  3.2  (continued) 
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W 

dG  = dx, 

Pc 


(3.2) 


where  Rs,  V,  and  dX  are  the  sheet  resistance  of  the  semiconductor 
layer  under  the  contact  hole,  the  width  of  the  contact  hole,  and  the 
length  of  a single  lump  of  the  distributed  model,  respectively.  The 
resistance  component  of  the  metal  layer  is  neglected  because  of  its 
small  resistivity.  The  current  flowing  to  the  contact  edge  in  the 
semiconductor  layer  will  travel  through  the  path  which  gives  the 
minimum  resistance. 

Muramann  and  Widmann  [6]  and  Berger  [7]  modeled  the  contact 
resistance  by  a Transmission  Line  Model  (TLM) , which  is  widely  used 
for  characterizing  the  horizontal-type  contact  [35-45].  To 
qualitatively  understand  the  basic  concept  of  the  contact  crowding, 
let  us  assume  that  d « V,  where  d is  the  length  of  a contact.  Then 

the  current  distribution  j(x)  across  the  contact  hole  can  be 
expressed  by 


lO  cosh/a(x  - d) 

j(x)  = 

V sinh/ad 


(3.3) 


where  IQ  is  the  total  current  through  the  contact  lead,  a = 

(Rs/pc)1/2,  and  Rs  is  the  sheet  resistance  of  the  semiconductor 
layer  under  the  contact. 

The  current  distribution  across  the  contact  hole  for  various 
values  of  y = x(Rs/pc)l/2  is  shovn  in  Fig,  3>3>  The  contact  with 

the  higher  Rs,  longer  d,  and  smaller  pc  has  heavier  contact 
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Fig.  3.3  Current  distribution  across  the  contact  hole  shown 
Fig.  3.2. 


m 
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crowding.  This  means  that  the  effective  resistance  of  this  kind  of 
contact  is  larger  than  the  value  from  the  simple  formula  Rc  = 
Pc/Area.  Note  also  that  for  y = 0 the  current  flows  uniformly  across 
the  contact.  This  situation  occurs  for  very  small  Rs  or  very  large 
PC- 


3-3  Model  for  the  Emitter  Resistance 
As  discussed  in  Section  3.2,  the  emitter  contact  of  a bipolar 
transistor  can  be  classified  into  one  of  two  categories  depending  on 
the  operating  condition.  For  the  emitter  contact  to  operate  as  an 
ideal  horizontal-type  contact,  the  emitter  current  should  flow  only 
into  the  one  side  of  the  contact.  This  cannot,  however,  occur  in  a 
real  transistor  as  is  explained  below.  Since  the  tendency  of  the 
emitter  current  to  flow  into  one  side  of  the  emitter  contact  is 
controlled  by  the  emitter  crowding,  we  have  to  accurately  analyze 
the  emitter  crowding.  Although  many  authors  treated  the  problem  of 
the  emitter  crowding  [8,  29],  their  approaches  are  not  applicable 
for  our  case  because  they  neglected  the  contribution  of  the  emitter 
resistance  to  the  emitter  crowding.  This  neglect  overestimates  the 
emitter  crowding,  because  the  negative  feedback  from  the  emitter 
resistance  is  not  taken  into  the  account. 

A single  transistor  can  be  modeled  as  a combination  of 
elementary  one-dimensional  transistors  [15].  By  adding  TLM  for  the 
emitter  contact  to  this  combination  of  elementary  transistors,  we 
obtain  the  equivalent  circuit  shown  in  Fig.  3.4.  For  a 
simplification  of  the  analysis,  we  make  the  following  assumptions: 
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(1)  the  current  gain  £ is  independent  of  position  and  current, 

(2)  the  base  conductivity  modulation  and  push-out  are 
negligible,  and 

(3)  the  emitter  has  a butted  structure. 

Although  the  above  assumptions  cause  some  error,  they  are  reasonable 
in  a sense  that  the  results  give  a correct  trend  for  the  emitter 
resistance  dependence  on  current,  pc,  sheet  resistances  of  the  base 
and  emitter,  and  geometry.  This  is  also  confirmed  experimentally. 
The  transistor  shown  in  Fig.  3.4  has  the  width  of  W and  the  length 
of  d.  The  elements  of  Fig.  3.4  are  defined  as  follows 

ARg  = Rgg.Ad/W 
Pc 

^EV  = 

Ad.W 

rSE  • Ad 

^EL  = 

V 

where  RSB  is  the  base  sheet  resistance,  RSE  is  the  emitter  sheet 
resistance,  Ad  = d/N,  and  N is  the  numbers  of  the  elementary 
transistors  used  in  modeling. 

3-4  Computer  Simulation 

In  the  previous  Section,  we  have  developed  an  equivalent 
circuit  model  for  the  purpose  of  modeling  the  emitter  resistance. 
Our  next  step  is  to  extract  the  effective  emitter  resistance  from 
the  model.  From  Fig.  3.4,  we  can  derive  a set  of  nonlinear 


(3.4) 

(3.5) 

(3.6) 
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Fig.  3.4  Equivalent  Circuit  of  sectioned  transistor  model. 

ARev  and  AREL  represent  the  resistances  from  interface 
layer  of  contact  hole  and  emitter  bulk(n+). 
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differential  equations  with  respect  to  position  x,  VBE(x),  and 
Ib(x)»  These  equations  can  be  solved  numerically.  It  is  easier, 
however,  to  analyze  the  circuit  of  Fig.  3.4  using  the  circuit 
simulator  SLICE,  which  is  a modified  version  of  SPICE  by  HARRIS. 
Another  advantage  of  using  SLICE  is  that  we  can  easily  sweep  all 
elements  of  the  circuit.  If  a transistor  is  divided  into  enough 
number  of  elementary  transistors,  the  results  from  SLICE  will  be 
equivalent  to  the  solution  of  the  differential  equations.  We 
confirmed  the  validity  of  using  SLICE  by  comparing  the  base 
resistance  values  for  the  single-base  contact  from  SLICE  with  the 
values  from  the  well-known  loe  current  formula  (RB  = RSBd/3W).  if  a 
transistor  is  divided  into  20  elementary  transistors,  the  deviation 
between  the  two  approaches  is  within  5%.  The  simulation  shown  below 
uses  20  elementary  transistors.  We  varied  four  basic  parameters,  and 
selected  the  values  of  the  parameters  as  close  as  possible  to  the 
values  for  the  real  devices.  The  ranges  for  four  parameters  are: 

RSB  = 3,000  to  11,000  Q/square 
Pc  = 10  to  200  8-pm2 
rSE  = 3 to  20  S/square 
d = 0.5  to  10  pm. 

The  results  from  the  SLICE  simulation  give  only  the  voltages  at 
the  nodes  and  the  currents  through  the  elements.  From  these,  we  have 
to  extract  the  resistance  values  for  the  emitter  and  base 
resistances.  The  power  dissipation  concept  [31]  was  used  in 
calculating  Rg  and  R£: 
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where  IB,  IE,  IB(i),  IEV(i),  and  IEL(i)  are  base  terminal  current, 
emitter  terminal  current,  the  base  current  of  an  elementary 
transistor  at  a certain  position  (i),  and  emitter  current  through 
the  interface  layer  of  the  contact  and  the  semiconductor  layer  in 
the  emitter  at  a certain  position  (i),  respectively.  The  vertical 
resistance  component  of  the  emitter  region  is  neglected  due  to 
heavily-doped  shallow  emitter.  The  first  term  of  (3.8)  represents 
the  emitter  resistance  component  from  the  interface  layer  of  the 
emitter  contact.  The  second  term  of  (3.8)  represents  the  emitter 
resistance  from  the  n+  layer. 


3*5  Results  and  Discussion 
3.5.1  Computer  Simulation  Results 

For  the  understanding  of  the  emitter  and  contact  crowding  of  a 
transistor,  we  first  consider  Fig.  3.5  that  shows  the  current 
distribution  for  VBE  = 0.85  V,  RSB  = 6000  52/square,  pc  = 30  52-pm2,  V 
= 10  urn,  and  d = 5 urn.  The  y-axis  of  this  graph  represents  the 
normalized  values  of  each  current  to  its  value  at  x = 0.25  pm.  The 
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Fig.  3.5  Current  distributions  for  the  transistor  in  Fig.  3.4. 

IEV»  IEL>  and  JB  are  the  currents  through  AREV,  AREL,  and 
ARb,  respectively.  Emitter  length  d = 5 pm,  specific 
contact  resistance  of  the  emitter  contact  pc  = 30  2-pm2, 
VfiE  = °*85  v»  the  base  sheet  resistance  RSB  = 6000 
2/square,  and  the  emitter  sheet  resistance  Rep  = 10 
2/square. 
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currents  at  x = 0.25  ym  correspond  to  the  first  transistor  in  Fig. 
3.4,  since  a transistor  with  d = 5 ym  is  divided  into  20  elementary 
transistors.  The  single  base  contact  is  located  at  x = 0 ym.  The 
base  current  IB(i)  decreases  with  increasing  x because  the 
elementary  transistors  far  away  from  the  base  contact  have  less 
forward  bias  at  the  base-emitter  junction  than  the  transistors  near 
the  base  contact  due  to  voltage  drop  in  the  base.  This  phenomena  is 
well-known  as  the  emitter  crowding.  As  the  emitter  current  crowds 
near  the  base  contact,  the  emitter  contact  starts  to  act  as  the 
horizontal-type  contact.  The  emitter  currents  IE(i)  Gf  the 

transistors  in  the  heavy  emitter-crowding  zone  will  not  pass 
entirely  through  AREV  connected  to  each  transistor.  If  the  entire 
emitter  current  of  each  transistor  flows  into  AREV  connected  to  it, 
the  emitter  resistance  becomes  large  due  to  the  significant 
contribution  from  AREV  = pc/(Ad  x V).  Thus,  IEV(i)  flows  more 
uniformly  compared  to  the  base  and  emitter  currents.  We  can  easily 
understand  this  mechanism,  if  we  check  the  distribution  of  the 
current  IEL  flowing  through  AREL.  The  emitter  current  of  Qj  in  Fig. 
3.4  passes  partially  through  AREV1,  the  rest  flows  through  the  next 
AREL.  This  tendency  continues  until  IE(i)  becomes  larger  than 
IEV(i).  Thus,  IEL  in  Fig.  3.5  increases  due  to  the  contribution  from 
the  rest  of  the  emitter  currents  of  the  previous  transistors  which 
have  enough  emitter  currents  in  supplying  IEV(i).  After  the  peak 
point,  IEL  decreases  to  support  insufficient  IE(i)  in  supplying 
IEV(i) • 
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A)  The  effect  of  the  emitter  crowding: 

Since  the  tendency  for  the  emitter  current  to  flow  into  one 
side  of  the  emitter  contact  depends  on  the  emitter  crowding,  the 
heavier  emitter  crowding  causes  heavier  contact  crowding.  The 
emitter  crowding  is  a function  of  RSB,  d,  and  VBE . Fig.  3.6  displays 
the  effective  emitter  resistance  (REFF),  normalized  to  low  current 
of  pc/Area,  as  a function  of  d.  REFF  increases  with  the 
emitter  length  and  becomes  larger  for  the  higher  RSB  at  the  same 
emitter  length.  The  dependence  of  REFF  on  VBE  is  depicted  in  Fig. 
3.7.  Figures  3.6  and  3.7  support  our  idea:  the  contact  crowding 

becomes  significant  for  heavy  emitter  crowding. 

B)  The  effect  of  the  contact  crowding: 

If  the  emitter  crowding  occurs,  we  have  to  also  consider  the 
contact  crowding  to  properly  model  REFF.  We  cannot  apply  the  theory 
developed  by  Murramann  and  Widmann  [6]  and  Berger  [7]  to  the  emitter 
contact  of  bipolar  transistor,  because  the  emitter  current  does  not 
flow  into  the  one  side  of  the  emitter  contact,  but  is  rather 
distributed  (see  Fig.  3.5).  It  is  hard  to  evaluate  the  emitter 
resistance  analytically.  Thus,  we  rely  on  computer  simulation.  To 

investigate  RFFF  as  a function  of  the  contact  crowding,  we  use  the 
following  parameters: 

pc  = 10  52-pm2  to  250  S-pm2  for  RSE  = 5,  10,15  52/square 
d = 0.5  pm  to  10  pin  for  pc  = 10,  60,  110  52-pm2. 
t,  we  consider  the  effect  of  pc  and  RgE  on  the  emitter 
resistance.  In  Fig.  3.8,  REFF/RE0  decreases  as  pc  increases.  For 
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Fig.  3.6  Effective  emitter  resistance  REFF  as  a function  of  emitter 
length.  Reff  is  normalized  to  the  low  current  emitter 
resistance  (oc/Area).  (VBE  = 0.9  V,  RSE  = 10  S/square, 

PC  = 50  S-pm2,  d = 10  pm) 
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Fig.  3.7  Effective  emitter  resistance  REFp  as  a function  of  VBE. 
(rSE  = *0  Q/square,  pc  = 50  C-pm*,  d = 10  pm, 
rSB  = 7000  S/square) 
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Fig.  3.8  Effective  emitter  resistance  as  a function  of  the  specific 
contact  resistance.  (VBE  = 0.9  V,  RSB  = 7000  S/square, 

RBe  = 10  Q/square,  d = 10  ym) 
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smaller  pc,  the  emitter  currents  of  the  transistors  in  the  current 
crowding  region  have  a strong  tendency  to  flow  through  AREV  instead 
of  traveling  further  across  the  emitter  region.  If  pc  is  large,  the 
emitter  currents  tend  to  spread  to  many  parallel  AREV  branches  to 
reduce  the  resistance  from  the  interfacial  layer.  REFF/RE0 
approaches  unity  at  high  pc  for  RSE  = 5 2/square,  but  for  RSE  = 10 
and  15  2/square  it  does  not  approach  unity  because  of  the  high 
lateral  resistance  component  from  the  emitter  region.  This  will  be 
discussed  in  the  next  paragraph.  Figure  3.9  shows  that  REFF/RE0 
increases  with  increasing  emitter  length  d.  This  is  because  for  long 
d we  have  both  heavy  emitter  and  contact  crowding. 

C)  The  components  of  the  emitter  resistance: 

In  Fig.  3.10,  we  illustrate  the  distribution  of  IEL  for  three 
various  values  of  VBE.  Since  the  emitter  current  flows  vertically  at 
low  current,  there  is  no  horizontal  potential  difference  in  the 
emitter  region.  At  low  VBE,  IEL  is  uniform  and  its  value  is 

negligible  compared  to  the  total  emitter  current.  For  example, 
IEL(Peak)  = 5.16  x 10-8  A and  IE  = 3.07  x 10-5  A for  VBE  = 0.75  V. 
Thus,  the  emitter  resistance  at  low  current  is  due  primarily  to  the 
contact  resistance.  But,  we  cannot  neglect  the  horizontal  component 
of  the  emitter  resistance  due  to  the  contact  and  emitter  crowding  at 
high  currents.  For  VBE  = 0.85  V the  peak  value  of  IEL  is  4.77  x 10-5 
A and  IE  = 9.73  x ltH  A.  Then,  the  contribution  of  AREL  is  not 
negligible  in  calculating  the  emitter  resistance  according  to  (3.8), 
because  the  sum  of  IEL(i)  is  comparable  to  IE.  The  two  components  of 
the  emitter  resistance  are  shown  in  Fig.  3.11.  It  is  clear  that  the 
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Fig.  3.9  Effective  emitter  resistance  as  a function  of  emitter 
length.  (Vgg  = 0.9  V,  Rgg  = 7000  Q/square,  RcE  = 10  2/ 
square,  pc  = 50  2-pm2) 
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POSITION  (MICROMETER) 


Fig.  3.10  Current  distribution  along  ARE£  as  a function  of  Vg£. 

(rSB  = 7000  2/ square,  RSE  = 10  2/square,  p = 50  Q-ym2, 
d = 10  ym) 
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Fig.  3.11  Variation  of  the  emitter  resistance  components  as  a 
function  of  VgE.  The  solid  line  is  for  the  contact 
resistance.  The  dotted  line  is  for  the  emitter  bulk 
resistance.  (Rgg  = 7000  Q/square,  RSE  = 10  Q/square, 
pc  = 50  Q-pm2,  d = 10  pm) 
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lateral  resistance  component  (AREL)  contributes  to  the  emitter 
resistance  at  high  VBE.  We  can  also  see  that  the  vertical  component 
(ARev)  increases  as  VBE  increases. 

3-5.2  Experimental  Results 

By  rearranging  equation  (2.6),  we  obtain 
AV  RB 


The  first  term  of  the  right  hand  side  of  (3.9)  either  decreases  [15] 
or  remains  constant  with  increasing  VBE  if  Rg/g  is  current 
independent  [47].  We  compared  two  extreme  cases  of  the  emitter 
contacts  represented  by  two  different  devices.  The  devices  are  a 
metalcontact  transistor  with  the  emitter  size  of  25  x 25  pm2  and  a 
polysilicon-contact  transistor  with  the  emitter  size  of  2 x 12  ym2. 
The  reason  for  choosing  these  devices  is  that  the  metalcontact 
device  can  have  a strong  emitter  and  contact  crowding  due  to  long  d 
and  low  pc  ( pc  of  the  metalcontact  device  is  smaller  than  that  for 
the  polysilicon  device).  On  the  other  hand,  the  polysilicon-contact 
device  has  a weak  emitter  and  contact  crowding  because  of  its  short 
d and  large  pc  ( pc  = 250  S2-um2  for  this  device  [27]).  Figure  3.12 
shows  the  experimental  results.  For  easy  comparison  of  the  two 
devices,  the  y-axis  displays  AV/IE  normalized  to  its  peak  value. 
AV/IE  for  the  metal-contact  device  continuously  increases  up  to  VBE 
= 1.0  V,  but  the  polysilicon-contact  device  has  the  opposite 
dependence.  The  AV/IE  increase  for  the  metalcontact  device  is  caused 
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Fig.  3.12  Experimental  results  for  two  transistors  with  different 

contact  materials  and  size.  The  emitter  contact  sizes  are 
2 x 12  pm  and  25  x 25  pm  for  polysilicon  and  metal 
contacts,  respectively.  For  easy  comparison,  the  y-axis 
(AV/Ie  = Re  + Rb/(1  + 3))  is  normalized  to  the  peak  value 
for  both  devices. 
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by  the  increase  of  RE  due  to  the  heavy  emitter  and  contact  crowding, 
which  agrees  with  the  simulation  results  in  Section  3.5.1.  The  AV/Iu 

ij 

decrease  after  VBE  = 1.0  V is  the  result  of  the  base  conductivity 
modulation  by  high  injection  phenomena.  Computer  simulation  results 
m Section  3.5.1  pointed  out  that  the  transistor  with  low  RSB  has 
only  small  small  increase  in  RE.  Thus,  the  metal-contact  transistor 
has  small  RE  increase  or  constant  RE  beyond  VBE  = 1.0  V. 

3.6  Conclusions 

Through  the  computer  simulation  based  on  our  new  model,  we  have 
found  that  the  emitter  resistance  increases  as  a function  of 
current.  The  dependence  of  the  emitter  resistance  on  physical 
parameters  was  also  investigated.  The  resistance  consists  of 
components  from  the  contact  hole  and  the  bulk  emitter  region.  The 
simulation  results  were  checked  by  measuring  devices  with  metal  and 
polysilicon  contacts.  The  experimental  results  show  that  the  emitter 
resistance  of  the  metal-contact  transistor  with  relatively  large 
emitter  length  increases  at  high  current.  This  conclusion  is 
important  for  the  design  of  analog  transistors  and  power 
transistors,  which  in  most  cases  have  long  emitter  with 


metalcontact . 


CHAPTER  FOUR 

MEASUREMENT  OF  THE  COLLECTOR  RESISTANCE 
4.1  Introduction 


The  delay  time  of  the  Emitter  Coupled  Logic  (ECL)  consists  of 
several  time-constants,  which  are  determined  by  the  series 

resistances  and  capacitances  of  a transistor  and  its  load 

resistance.  The  delay  time  can  be  expressed  as  a the  sum  of  eight 
time-constants  [1].  The  resistance,  which  includes  the  collector 
series  resistance  (Rc)  and  the  actual  load  resistance,  contributes 
to  five  out  of  eight  time-constants.  Thus,  Rc  plays  an  important 
role  in  limiting  the  speed  of  digital  circuits.  Another  important 
aspect  of  RC  in  the  digital  applications  is  that,  together  with  the 
emitter  resistance  and  the  internal  saturation  voltage  VCE(sat),  it 
determines  the  value  of  the  logical  low  state.  The  transition 
frequency  (fT)  of  a transistor  is  used  as  a figure  of  merit  in 
estimating  the  frequency  response  of  a transistor  in  the  analog 
application.  The  fT  is  also  a strong  function  of  the  Rc,  together 
with  other  parameters  of  a transistor. 

Therefore,  Rc  is  a crucial  parameter  in  estimating  the 
switching  speed  and  frequency  response  of  the  bipolar  transistor.  It 
is  also  one  of  the  most  difficult  parameters  to  measure  due  to  its 
strong  dependence  on  bias  and  its  distributed  nature. 

The  collector  resistance  of  an  integrated  bipolar  n-p-n 
transistor,  shown  in  Fig.  4.1,  consists  of  the  vertical  portion  Rcv 
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Fig.  4.1  Device  cross-section  showing  the  different 
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of  the  n type  epitaxial  region  between  the  edge  of  the  base- 
collector  space-charge  region  and  the  n+  subcollector,  the  lateral 
portion  of  the  n+  subcollector,  and  another  vertical  portion  of  the 

n+  plug  between  the  subcollector  and  the  collector  contact.  The 
contact  resistance  has  to  be  also  considered. 

The  heavily  doped  portions  of  Rc  plus  the  contact  resistance 
are  current  independent  and  will  be  lumped  together  into  one 
resistance  RCL.  The  value  of  the  vertical  portion  Rcv  - 
Pepitepi/ABE*  ^ a complicated  function  of  current,  voltage,  and 
doping.  At  low  currents,  the  thickness  of  quasi-neutral  epitaxial 
collector  tepi  depends  on  the  base-collector  junction  voltage,  that 
determines  the  junction  space-charge  region  width.  At  high  currents, 
the  base-collector  junction  may  become  forward-biased  resulting  in  a 
conductivity  modulation  of  the  epitaxial  region  that  lowers  Rcv.  The 
epitaxial  layer  will  be  also  modulated  by  base-pushout  at  high 
currents.  To  accurately  model  Rc,  ve  have  to  measure  it  as  a 
function  of  current  and  biasing  conditions.  We  present  here  two  DC 
methods  for  measuring  RCL  and  Rcv.  The  methods  are  demonstrated  on 
advanced  MOSAIC  III  self-aligned  polysilicon  transistors. 

The  existing  methods  for  measuring  the  Rc  will  be  reviewed  in 
Section  4.2.  The  basic  concept  of  our  new  method  for  measuring  Rc 
will  be  discussed  in  Section  4.3.  In  Section  4.4,  the  detailed 
measurement  procedure  will  be  introduced  and  the  results  for  various 
samples  are  shown  and  discussed. 
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4*2  Review  of  the  Conventional  Methods 
A short  critical  review  of  the  published  methods  to  measure  Rc 
is  presented  here.  Detailed  description  is  presented  elsewhere  [24]. 

The  high-current  value  Rc(sat)  and  the  low-current  value 
Rc(normal)  of  Rc  can  be  estimated  from  the  IC-VCE  characteristics  of 
a transistor.  Rc(normal)  can  be  estimated  from  the  slope  of  a line 
drawn  through  the  "knees"  of  the  IC-VCE  curves  for  different  values 
of  the  base  current  IB.  It  is  very  difficult,  however,  to  determine 
the  "knees"  for  each  curve.  Also,  Rc(normal)  from  this  measurement 
cannot  represent  the  values  for  the  entire  forward-active  region, 
since  Rc(normal)  depends  strongly  on  the  collector  current  and  the 
collector-base  reverse  bias  due  to  the  space-charge  region 
modulation  in  the  lightly  doped  collector.  To  obtain  Rc(sat),  we 
analyze  the  expression  for  VCE(sat)  given  by 


VCE(sat) 


kT  n + (IC/IB)(1  - Or)' 

In  — 

q “R  (1  - (Ic2  > 


+ IrRe  + IcRC» 


(4.1) 


where  or  is  the  common-base  reverse  current  gain,  fjp  is  the  common- 
emitter  forward  current  gain,  and  the  other  parameters  in  (4.1)  have 
their  usual  meanings.  From  (4.1),  we  can  obtain  RE  + Rc  by  plotting 
VCE(sat)  as  a function  of  Ic,  provided  that  the  first  term  remains 
current  independent.  The  effect  of  the  first  term  of  (4.1)  can  be 
reduced  using  the  forced  beta  method,  in  this  method,  we  measure  the 
IB-VCE  plots  for  the  several  values  of  Ic.  Ve  then  find  AVCE  as  the 
difference  between  two  points  on  two  different  curves  for  the  same 
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ratio  0 - IC/IB.  If  we  further  assume  that  and  gF  are  independent 
of  Ic,  and  that  RE  « Rc,  then  from  (4.1)  Rc  = AVCE/AIC. 

A method  which  uses  the  parasitic  p-n-p  transistor  of  the 
integrated  n-p-n  transistor  to  measure  Rc  was  proposed  [48].  At  high 
currents,  the  collector  current  induces  forward  bias  at  the 
collector-base  junction  of  the  n-p-n  transistor  due  to  the  voltage 
drop  on  the  collector  series  resistance.  The  forward-biased 
collector-base  junction  of  the  n-p-n  transistor  behaves  as  the 
emitter  of  the  p-n-p  transistor.  The  holes  injected  into  the  epi 
that  cross  the  subcollector  are  collected  by  the  substrate  and 
constitute  a substrate  current,  Is. 

Since  the  base  of  the  p-n-p  is  distributed  across  Rcv  (Fig. 
4.1),  the  question  is  where  to  attach  the  base  terminal  for  the 
purpose  of  calculating  the  base-collector  potential  that  controls 
IS*  Based  on  numerically  calculated  field  profile  in  the  epi,  the 
base  terminal  of  the  p-n-p  transistor  was  positioned  at  the  epi- 
subcollector  interface  in  the  reference  [48].  As  will  be  shown 
later,  the  creation  of  the  base  terminal  like  this  causes  an  error 
in  determining  of  Rc  from  Is.  We  also  show  that  Ig  is  controlled 
directly  by  the  physical  base-collector  junction  of  the  n-p-n 
transistor. 

For  accurate  modeling  of  a bipolar  transistor,  the  two 
components  of  Rc  (Rcv  and  Rcv)  have  to  be  separated  and  their  bias 
dependence  has  to  be  determined.  Neither  of  the  above  methods  can  do 
that.  Our  new  method  clearly  separates  the  two  components,  and  also 
yields  their  bias  dependence.  In  addition,  a new  procedure  for 
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separating  RE  from  Rc  using  the  forced  beta  method  will  be 
presented. 


*♦  • ->  myaxcax  noaei. 

A careful  Investigation  of  the  cross-section  of  the  integrated 
n-p-n  transistor  reveals  the  presence  of  a parasitic  p-n-p 
transistor,  as  shovn  in  Figs.  4.2  (a)  and  4.3  (a).  In  these  figures, 
RCL  represents  the  resistance  from  the  collector  contact  to  the 
collector-substrate  junction,  and  Rcv  is  resistance  of  the  quasi- 
neutral epitaxial  layer.  The  parasitic  p-„-p  transistor  can  be 
turned  on  by  the  voltage  drop  on  the  collector  resistance  of  the  n- 
p-n  transistor.  The  emitter  of  the  p-n-p  transistor  can  be  either 
the  base  or  the  substrate  of  the  n-p-n  transistor,  depending  on  bias 
conditions.  Be  define  the  upward  operation  of  the  p-n-p  transistor, 
when  the  substrate-collector  junction  of  the  n-p-n  transistor  is 
forward  biased,  and  the  downward  operation  when  the  base-collector 
junction  of  the  n-p-n  transistor  is  forward  biased. 


4‘3,1  HEvard  Operation  of  the  p-n-p  Transistor 

If  the  reverse  bias  at  the  collector-base  junction  is  larger 
than  the  voltage  drop  on  the  RCL  and  Rcv,  the  collector-base 
junction  cannot  turn  on.  Under  this  condition,  however,  the 
substrate-collector  junction  can  turn  on  due  to  the  voltage  drop  on 
Rci-  The  substrate  operates  then  as  the  emitter  of  the  p-n-p 
transistor.  Since  the  carrier  injection  into  the  base  of  the  p-n-p 
P-n-p  transistor  occurs  at  the  substrate-collector  junction  in  this 
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EMITTER 


COLLECTOR 


(a) 

Fig.  4.2  Illustration  of  upward  operation  of  the  parasitic  p-n-p 
transistor. 

(a)  Device  cross-section  for  upward  operation  of  the 
parasitic  p-n-p  transistor.  is  the  resistance  of 
the  n epi-layer,  RqV  is  the  resistance  of  the  n+ 
subcollector  and  plug. 

(b)  Equivalent  circuit  for  upward  operation  of  the 
parasitic  p-n-p  transistor. 
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Fig.  4.2  (continued) 
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operation,  the  emitter  current  of  the  p-n-p  transistor  is  controlled 
by  the  effective  emitter-base  junction  voltage  (VEB)  of  the  p-n-p 
transistor  determined  by  the  substrate  and  collector  terminal 
voltage  minus  the  voltage  drop  on  RCL  (see  Fig.  4.2  (a)). 

For  detailed  analysis,  we  show  the  equivalent  circuit  for  the 
upward  operation  of  the  p-n-p  transistor  in  Fig.  4.2  (b).  We  assume 
that  the  RCL  is  independent  on  the  current  levels,  because  the 
buried  layer  is  usually  heavily-doped.  The  voltage  drop  at  the  Rs  is 
negligible  in  comparison  with  other  components  if  the  R^  is 
obtained  at  the  relatively  low  substrate  current,  for  example,  Is  < 

1 PA.  Then,  VEB  of  p-n-p  for  two  different  substrate  voltage  (Vs), 
is  expressed  as 

VEB1  = VS1  - IS1RS  + IC1RCL  _ vc  (4.2) 

VEB2  = VS2  ~ iS2rS  + iC2rCL  ~ VC*  (4.3) 

For  the  same  Is,  VEB  will  be  the  same.  Thus,  equaling  (4.2)  and 
(4.3)  gives 


rCL 


VS2  - VS1 
IC1  " ZC2 


(4.4) 


The  rCL  component  of  Rc  is  found  uniquely  from  (4.4).  This  is  a 
great  advantage  in  comparison  with  other  methods  where  different 
components  of  the  series  resistances  are  always  coupled  together. 
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4-3.2  Dovnvard  Operation  of  the  p-n-p  Transistor 

For  the  downward  operation  of  the  p-n-p  transistor,  the  base- 
collector  junction  is  forward-biased,  while  the  substrate-collector 
junction  is  reverse-biased.  Then,  the  potential  right  below  the 
collector-base  junction  controls  the  emitter  current  of  the  p-n-p 
transistor,  because  the  injection  occurs  at  this  junction.  Thus,  the 
base,  collector,  and  substrate  of  the  n-p-n  transistor  will  behave 
as  the  emitter,  base,  and  collector  of  the  p-n-p  transistor,  as 
shown  in  Figs.  4.3  (a)  and  4.3  (b).  The  holes  injected  from  the  p+ 
base  transverse  through  the  epitaxial  layer  and  the  subcollector  and 
constitute  the  collector  current  of  the  p-n-p  transistor.  This 
current  is  also  the  substrate  current  of  the  n-p-n  transistor  (see 
Fig.  4.3  (b)).  The  collector  current  of  the  p-n-p  transistor  is 
controlled  by  the  potential  VEB  at  top  p-n  junction  determined  by 

the  terminal  voltages  VB  and  Vq  minus  the  voltage  drop  on  both 
and 

VEB1  = VB  ~ rB1iB1  + RC1IC1  ~ VC1  (4.5) 

VEB2  = VB  ~ rB2iB2  + rC2iC2  ~ vC2»  (4.6) 

where  the  subscripts  1 and  2 correspond  to  VC1  and  vc2, 
respectively.  The  collector  current  of  a transistor  remains  fixed 
when  VEB  and  VCB  are  fixed.  Thus,  for  the  same  Ig,  VEB  in  (4.5)  and 
(4.6)  will  be  the  same: 

VC1  " VC2  AVC 

IC1  - XC2  AIC 


RC  = RCV  + rCL  ~ 


(4.7) 
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EMITTER 


(a) 

Fig.  4.3  Illustration  of  the  downward  operation  of  the  parasitic 
p-n-p  transistor. 

(a)  Cross-section  for  downward  operation  of  the  parasitic 
p-n-p  transistor. 

(b)  Equivalent  circuit  for  the  downward  operation  of  the 
parasitic  p-n-p  transistor. 
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Fig.  4.3  (continued) 
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Two  assumptions  were  used  in  (4.7).  The  first  assumption  is 
that  the  voltage  drop  on  the  base  resistance  for  two  different 
values  of  Vc  is  almost  the  same.  This  assumption  can  be  verified 
experimentally.  For  example,  the  voltage  difference  (IbiRb1  - 
iB2rB2)  for  the  transistor  with  the  emitter  size  of  2 x 12  pm2  is 
only  about  2 mV  for  AVq  = 30  mV.  The  second  assumption  is  that  Rq  is 
constant  for  small  variation  of  Vc  (AVC  < 50  mV).  The  small 
variation  of  Vc  may  modulate  the  base  width  of  the  p-n-p  transistor, 
but  this  modulation  is  negligible  for  small  AVq. 

In  one  previous  work  [48], the  potential  controlling  Is  was 
chosen  to  be  at  the  epi-subcollector  junction.  This  corresponds  to 
replacing  Rc  in  (4.5)  and  (4.6)  by  RCL.  This  neglects  the  voltage 
drop  on  Rcv.  We  show  in  the  next  section,  that  the  value  obtained 
for  the  downward  operation  mode  includes  Rcv. 

4.4  Experimental  Results 
4.4.1  Measurement  of  the  Collector  Resistance 

The  measurements  were  performed  on  advanced  self-aligned 
polysilicon-contact  transistors.  Figure  4.4  shows  the  current- 
voltage  characteristics  for  transistor  with  the  emitter  mask  size  of 
1.75  x 4.75  pm2.  Figure  4.5  shows  the  Ic  and  Is  versus  VBE  curves 
for  this  transistor.  The  measurements  were  done  for  fixed  VB  = 0 V 
and  Vc  = 1 V,  while  Vs  was  stepped  by  A Vs  = 30  mV  from  Vs  = 1 V.  The 
polarity  of  Is  is  positive  (Is  flows  into  the  substrate),  which 
confirms  that  the  bottom  p-n  junction  is  forward-biased.  Using  the 
IS  - VBE  curves,  we  find,  for  the  same  value  of  Is,  AIC 
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Fig.  4.4  Gummel  plots  for  the  transistor  with  the  emitter  mask  size 
ABE  = 1.75  x 4.75  ym^.  (Vg^  = 0 V) 
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corresponding  to  AVS  shown  in  Fig.  4.5.  Then,  RCL  = AVS/AIC  from 
(4.4).  Figure  4.6  shows  the  measured  values  for  RCL  as  a function  of 
^BE  ^or  three  transistors  with  the  emitter  mask  sizes  of  1.75  x 4.75 
pm  , 1.75  x 8 ym  , and  3x8  pm2,  respectively.  As  expected,  RCL  is 
almost  constant  in  the  range  from  VBE  = 0.9  V to  1.1  V.  We  note, 
that  Rcl  for  1.75  x 8 pm2  is  larger  than  that  for  3x8  pm2.  This 
can  be  explained  in  the  following  way.  The  collector  current  flows 
laterally  through  the  buried  layer  in  the  extrinsic  region  of  the 
transistor.  But,  the  current  flows  mostly  vertically  in  the  buried 
layer  under  the  intrinsic  transistor  area.  The  transistor  with  the 
larger  emitter  size  will  have  the  smaller  contribution  to  RqL  from 
the  intrinsic  buried  layer  in  comparison  with  the  smaller  emitter 
transistor,  because  the  resistance  is  inversely  proportional  to  the 
area  perpendicular  to  the  current  flow.  Thus,  even  if  the  spacing 
between  the  the  collector  contact  and  the  intrinsic  region  of  the 
transistors  is  the  same,  they  can  have  different  RCL  due  to  the 
reason  mentioned  above. 

The  measured  curves  for  the  downward  operation  of  the  p-n-p 
transistor  are  shown  in  Fig.  4.7.  The  collector-substrate  junction 
is  reverse-biased  with  Vs  « - 1 V,  and  VB  = 0 V.  We  measure  Ic  and 
Ig  as  a function  of  VBE  for  two  different  values  of  Vq  = - 0.285  V 
and  - 0.315  V.  Because  Ic  is  a strong  function  of  the  base  width, 
the  base-width  modulation  for  AVq  should  be  considered.  The  base  of 
the  parasitic  p-n-p  consists  of  the  epitaxial  layer  and  the  buried 
layer  with  total  base-width  of  about  0.6  pm.  Thus,  the  base-width 
modulation  for  small  AVC  = 30  mV  is  negligible.  This  fact  is 
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4.5 


Ic  and  Is  as  a function  of  Vgg  for  transistor  with 
ABE  = 1.75  x 4.75  pm^.  The  bottom  Ig  is  for  Vs  = 1 V, 
upper  Ig  curve  is  for  Vg  = 1.03  V. 
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4.6  Lateral  collector  series  resistances  for  three  different 
transistors  as  a function  of  VBE.  ABE  represents  the 
emitter  mask  size. 
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Fig.  4.7  Ic  and  Is  for  the  transistor  with  ABE  = 1-75  x 4.75  ym2 

versus  VBE.  Measurement  conditions  are  VB  = 0 V,  Vg  = -1  V, 
Vq  = -0.285  V for  the  top  Ig  curve,  and  Vq  = -0.315  V for 
bottom  Ig  curve. 
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confirmed  by  the  fact  that  the  Ic  - VBE  curves  in  Fig.  4.7  for  both 
values  of  Vc  are  nearly  identical.  The  total  collector  resistance  Rc 
= RCV  + rCL  ^ obtained  from  Fig.  4.7  in  a similar  way  as  for  the 
previous  case  of  the  upward  operation  by  measuring  AIC  corresponding 
to  AVC,  and  then  using  (4.7). 

To  investigate  the  Rc  dependence  for  VBE  <0.9  V,  the  current 
offset  ( = 1 pA)  can  be  subtracted  from  the  Is  values  measured  for 
higher  biases.  It  is  not  easy  to  measure  Is  in  the  pA  range  due  to 
the  leakage  current  between  the  leads  of  the  package.  Cleaning  the 
package  and  the  plug-in  of  the  shielded  box  of  the  HP4145  analyzer 
with  alcohol  dramatically  reduces  the  leakage  currents.  To  obtain 

the  flat  offset  currents  of  Is,  the  hold-time  of  HP4145  was  set  to 
20  seconds. 

The  measured  dependencies  of  Rc  on  VBE  and  Ic  for  two 
transistors  are  shown  in  Figs.  4.8  (a)  and  (b)  and  in  Fig.  4.9.  At 
high  currents,  Rc  saturates  at  Rc  = RCL,  which  agrees  very  well  with 
the  Rcl  dependence  found  in  Fig.  4.6.  For  smaller  currents,  Rc 
sharply  increases  because  of  the  contribution  from  the  epitaxial 
layer  of  Rcv.  The  very  sharp  transition  from  the  low-current  value 

to  the  high-current  value  is  caused  by  the  conductivity  modulation 
of  the  epitaxial  layer. 

To  investigate  the  dependence  of  Rc  on  bias  in  more  physical 
detail,  we  performed  numerical  simulations  using  MEDUSA  simulator. 
Figs.  4.10  (a)  and  (b)  show  the  simulated  dependence  of  Rc  on  VBE 
and  Ic  for  Vc  = - 0.3  V for  a transistor  with  1.75  x 4.75  pm2 
emitter  size.  The  doping  profiles  and  geometry  for  the  simulated 
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Fig.  4.8  Total  collector  resistances  versus  (a)  VBE  and  (b)  Ic 
for  three  different  values  of  Vc.  ABE  = 1.75  x 4.75 
pm2,  Vs  = -1  V. 
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Fig.  4.8  (continued) 
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Fig. 


4.9  Total  collector  resistance  versus 
values  of  Vc.  ABE  = 1.75  x 8 ym2, 


VgE  for  two  different 
Vs  = -1  V,  1.5  V. 
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transistor  are  similar  to  that  measured  in  Fig.  4.8.  The  qualitative 
agreement  between  Figs.  4.8  and  4.10  is  excellent. 

Rc  in  Fig.  4.10  is  calculated  from  the  quasi-Fermi  potential 
difference  between  the  collector  contact  and  the  edge  of  the  base- 
collector  space-charge  region.  The  boundary  between  the  quasi- 
neutral collector  and  the  base-collector  space-charge  region  is 
rather  easy  to  determine  at  small  and  moderate  biases  from  the 
electric  field  profile.  This  gets  more  difficult  at  biases  where  the 
quasi-neutral  base  pushes  out  into  the  epi  region  and  the  electron 
concentration  is  comparable  or  larger  than  the  collector  donor 
concentration.  It  is  easier  to  to  approximate  the  boundary  by  using 
the  difference  between  the  electron  and  doping  concentrations  and 
identify  the  boundary  at  a position  at  which  n = ND.  This  criterion 
was  used  to  calculate  Rq  displayed  in  Fig.  4.10. 

Figure  4.11  shows  the  detailed  simulated  dependencies  of  the 
excess  electron  concentration,  electric  field,  and  Rc  on  the 
distance  from  the  emitter  surface  for  VBE  = 0.87  V (Ic  = 0.82  mA). 
Rc  is  measured  from  the  collector  contact.  Also  shown  is  the  log  n - 
log  ND  curve.  Note  that  n = ND  for  x = 0.52  pm  where  Rc  = 203  52.  For 
small  currents,  Rc  is  almost  constant  and  equal  to  Rcv.  At  VBE  = 
0.87  V,  the  injected  electron  concentration  approaches  the  doping  in 
the  epi  layer  (see  Fig.  4.11).  As  a result,  Rcv  (and  Rc)  rapidly 
drops,  as  shown  in  Fig.  4.10  (a).  Simulation  shows  that  for  VBE  > 
0.89  V,  n > Nd,  and,  hence,  Rc  drops  to  a low  value  of  Rc  = RCL. 
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(a) 

Fig.  4.10  Simulated  total  collector  resistance  from  MEDUSA  versus 
(a)  VBE  and  (b)  Ip  for  the  transistor  with  Aop  = 

1.75  x 4.75  urn2. 
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Fig.  4.10  (continue) 
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Fig.  4.11  Simulation  results  from  MEDUSA  for  Vg£  = 0.87  V, 
^BC  = V,  Agg  = 1.75  x 4.75  ym^. 
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The  sharp  transition  between  the  low-current  and  high-current 
values  of  Rc  shown  in  Figs.  4.8  and  4.9  strongly  depends  on  the 
external  collector  bias  Vc.  The  internal  base-collector  potential  is 
given  by  the  external  voltages  Vc  and  VB  minus  the  voltage  drop  on 
the  collector  resistance.  As  Vc  in  Fig.  4.8  (a)  changes  from,  for 
example,  0 V to  + 0.3  V,  the  base-collector  junction  becomes  more 
reverse-biased.  Thus,  it  takes  larger  Ic  to  induces  enough  voltage 
drop  RCIC  to  forward  bias  the  junction.  The  critical  collector 
current  for  the  onset  of  the  base  push-out  also  increases  with  the 
increasing  base-collector  reverse  bias  in  ohmic  saturation  and  quasi 
saturation  [11-12].  The  transition  between  the  low-  and  high-Rc  is, 
thus,  moved  to  higher  currents  as  Vc  increases  (Fig.  4.8(b)).  Since 
the  collector  resistance  after  the  onset  of  the  base  push-out  is  a 
strong  function  of  the  high  injection,  it  is  not  easy  to  clearly 
explain  the  values  for  the  plateau  in  Fig.  4,8. 

4*4*2  Measurement  of  the  Emitter  Resistance 

As  discussed  in  Section  4.2,  the  slope  of  the  Ic  - VCE  curve  in 
saturation  is  proportional  to  the  sum  RCL  + re.  As  the  device  size 
gets  smaller,  RE  becomes  larger  because  RE  is  inversely  proportional 
to  the  emitter  area,  also  becomes  smaller  due  to  the  shorter 
distance  from  the  collector  contact  to  the  intrinsic  collector-base 
junction.  Nevertheless,  RE  may  become  comparable  to  RCL.  A very 
accurate  value  for  RCL  can  be  obtained  from  Is  of  the  upward 
operation  mode  of  the  p-n-p  transistor,  as  was  described  above. 
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Below  we  describe  a new  method  for  extraction  of  RE  based  on 
measurement  of  RCL  combined  with  the  foced  beta  method. 

Before  describing  the  detailed  procedure,  we  point  out  the 
problem  encountered  during  the  measurement  of  Rc(sat)  + (1  + l/g)RE 
by  taking  the  slope  of  the  IC-VCE  plot.  Ve  found  that  the  slope 
varies  in  a wide  range.  For  example,  for  the  transistor  with  the 
emitter  size  of  1.75  x 4.75  ym2  the  slope  varies  in  between  68  Q and 
81  S,  depending  on  the  points  at  which  the  slopes  are  taken.  The 
reason  for  the  various  slopes  is  that  the  internal  saturation 

voltage  VCE(sat)  is  not  negligible  in  comparison  with  Rc(sat)  + (1  + 
1/0)Re. 

Thus,  we  chose  the  forced  beta  method  that  overcomes  the 
problem  with  the  variable  slope.  For  IC/IB  fixed,  we  derive  from 
(4.1) 

AVCE  AIb 

— = Rc(sat)  + (1  + )RE.  (4. H) 

AIC 

The  SUm  of  the  collector  and  emitter  resistances  can  be  obtained 
accurately  from  (4.11).  This  method  is  usually  used  to  measure 
Rc(sat)  by  either  neglecting  RE  or  estimating  RE  by  the  open- 
collector  method  [24].  It  was  shown  [27],  however,  that  the  open- 
collector  method  is  not  reliable  and  grossly  overestimates  RE.  We 
can  use  (4.11)  to  extract  RE,  provided  that  Rc(sat)  can  be  measured 
independently  and  accurately.  Under  the  heavily  saturated  operating 
conditions  Rc(sat)  = RCL,  which  can  be  independently  measured  from 
the  upward  mode  of  the 


p-n-p  transistor. 
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4.12  Ig  versus  Vq  with  Iq  as  a parameter  for  transistor  with 
Abe  = 1*75  x 4.75  pm^.  The  collector  currents  are  from 
left  to  right:  0 mA,  0.2  mA,  0.4  mA,  and  0.6  mA. 
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In  Fig.  4.12,  we  illustrate  a typical  IB  -VCE(sat)  plot  for  the 
transistor  with  the  emitter  mask  size  of  1.75  x 4.75  ym2.  We 
investigated  the  variation  of  AVCE/AIC  as  function  of  current.  The 
results  are  summarized  in  Table  4.1  calculated  from  (4.11)  and  using 
the  Rcl  data  from  Fig.  4.6  together  with  the  values  for  Rg.  The 
variation  of  Rg  is  within  10  %.  Although  RCL  is  about  six  times 
larger  than  Rg,  the  two  terms  in  (4.11)  are  about  the  same  since 
MB/AIC  ~ 4-6  (see  Table  4.11).  The  accuracy  of  the  extracted  value 
for  Rg  depends  on  the  accuracy  of  RCL.  The  new  method  presented  here 
is  sufficiently  accurate  for  this  purpose.  The  values  for  Rg  for 
three  transistors  are  summarized  in  Table  4.2.  By  subtracting  0.4  ym 
wide  sidewall  from  the  emitter  mask  size,  we  obtain  the  effective 
emitter  area.  From  the  effective  emitter  area  and  the  values  for  Rg 

from  Table  4.2,  we  obtain  the  average  value  for  the  specific  contact 
resistivity  of  pc  = 50  B-ym2. 

4.5  Conclusions 

New  DC  methods  to  measure  the  collector  resistance  Rc  and 
emitter  resistance  Rg  of  integrated-circuit  bipolar  transistors  are 
is  presented.  The  method  for  Rc  is  based  on  monitoring  the  substrate 
current  of  the  parasitic  p-n-p  transistor  linked  with  the  n-p-n 
intrinsic  transistor.  The  p-n-p  transistor  is  operated  with  either 
the  bottom  or  top  p-n  junction  forward-biased.  This  allows  to 
separate  the  two  components  of  Rc.  Rg  is  obtained  from  the  measured 
lateral  portion  of  Rc  and  collector-emitter  saturation  voltage. 
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Table  4.1  Dependence  of  RE  on  Ic  for  the  transistor  with 
the  emitter  mask  size  of  1.75  x 4.75  pm2. 

RE  was  obtained  from  (4.8)  using  RCL  = 64  2 from 


Ic 

IB 

w 

o 

> 

avce/ aic 

rE 

(mA) 

(mA) 

(mV) 

(2) 

(2) 

0.2 

0.9 

38 

0.3 

1.35 

51 

130 

12 

0.2 

0.8 

36 

0.4 

1.6 

60 

120 

11.2 

0.2 

1 

38 

0.4 

2 

64 

130 

11 

0.2 

1.2 

41 

0.4 

2.4 

68 

135 

10.14 
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Table  4.2  Summary  of  RE  for  three  emitter  mask  sizes. 


abe 

(um2) 

rCL 

(2) 

rE 

(2) 

1.75  x 4.75 

64 

11  - 12 

1.75  x 8 

41 

7 ~ 7.58 

3x8 

37 

4 ~ 4.5 
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Examples  of  measurements  on  advanced  self-aligned  transistors  with 
polysilicon  contacts  are  shown  and  compared  with  numerical 
simulations.  The  results  show  very  strong  dependence  of  Rc  on  the 
base-emitter  and  base-collector  voltages  of  the  n-p-n  transistor. 
The  bias  dependence  of  Rc  is  due  to  the  conductivity  modulation  of 
the  epitaxial  collector.  From  the  measured  RE  we  obtain  a value  for 

the  specific  contact  resistance  for  the  polysilicon  emitter  contact 
of  pc  =50  Q-ym2. 


CHAPTER  FIVE 

MEASUREMENTS  OF  THE  AUGER  COEFFICIENTS 
5.1  Introduction 


There  are  large  discrepancies  of  the  Ambipolar  Auger 

coefficients  CA  = Cn  + Cp  from  one  investigator  to  another.  The 

values  from  Dziewior  and  Schmid  [49]  and  Svantesson  and  Nilsson  [50] 

for  the  plasma  density  of  above  1 x 1018  cnr3  (3.8  x 10“31  - 3.4  x 
31  6 1 

10  cm  s-  ) are  almost  same.  Dziewior  and  Schmid  extracted  the 
carrier  lifetime  from  the  decay  time  of  the  band-to-band 

recombination  radiation  after  pulsed  excitation  by  a mode-locked 
cavity-dumped  Ar+-ion  laser.  Svantesson  and  Nilsson  used  almost  the 
same  method  except  for  measuring  the  recombination  radiation  in  a 
direction  perpendicular  to  the  laser  beam  ( YAG : Nd ) . The  absorption 
of  infrared  probe  radiation  by  nonequilibrium  carriers,  excited  with 
light  from  a Neodymium  laser,  was  used  to  measure  CA  [51].  The  value 
for  CA  from  this  experiment  is  1.5  x 10"30  cm^"1  for  the  plasma 
density  above  5 x 101?  cnr3.  The  CA  from  the  photo-conductivity  and 
the  dynamic  grating  experiment  is  1.1  x lO'30  cmV1  for  the  plasma 
density  below  1 x lO*8  cnr3  [52].  In  the  experiments  based  on 
measuring  the  lifetime,  the  surface  recombination  of  the  sample  is 
always  a problem. 

Two  approaches  to  solve  this  problem,  by  a chemical  preparation 
method  [53]  or  by  growing  good  quality  oxide  layer  on  the  surface  of 
the  device  [54]  were  reported.  A strobe  lamp  excites  carriers  in  a 
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semiconductor  wafer  immersed  in  hydrofluoric  acid,  and  their  density 
decay  is  then  monitored  by  an  inductively  coupled  rf  bridge  [53]. 
The  value  from  this  measurement  is  2 x 10~30  cm6s"1  for  the  plasma 
density  in  the  range  of  3 x lO1^  to  1017  cm"3.  By  monitoring  the 
open-circuit  voltage  of  a highly  injected  solar  cell  under  steady- 
state  and  transient  conditions,  the  CA  was  extracted  [54].  The  value 
reported  from  this  experiment  is  1.66  x 10~30  cmV1  for  the  plasma 
densities  in  the  range  of  1 x 1016  to  1.2  x 1017  cm"3. 

In  this  chapter,  we  present  measurements  of  Auger  lifetime  in 
the  carrier  density  range  from  1017  to  5 x 1017  cm"3.  Our 

experimental  method  is  fundamentally  different  from  other  methods 
[49-54],  since  we  measure  the  electroluminescence  (EL)  as  a function 
of  the  current  through  a n+/p-/p+  diode.  To  extract  the  Auger 
recombination  coefficient,  the  experimental  results  are  compared 
with  numerical  calculations  that  self-consistently  solve  the 
transport  equations  across  the  p~  base  of  the  diode.  The  connection 
between  experimental  and  numerical  results  is  made  possible  by 
normalizing  the  electroluminescence  signal  at  high  currents  with 
respect  to  the  same  signal  under  low  injection.  For  carrier  plasma 
densities  between  lO1*  and  1017  cm"3,  recombination  in  the  n+  and  p+ 
layer  dominates  the  current  and  the  relation  between 
electroluminescence  and  current  is  linear.  In  this  density  range, 
fitting  experimental  results  with  numerical  calculations  reveals  the 
saturation  current  density  of  the  heavily  doped  regions.  For  plasma 
densities  between  iO37  and  5 x 1037  cnr3,  Auger  recombination 
reduces  the  base  ambipolar  lifetime  well  below  its  Shockley-Read- 


99 


Hall  value.  The  base  recombination  becomes  comparable  to  the 
recombination  in  the  heavily  doped  regions  and  the 

electroluminescence  intensity  starts  increasing  sublinearly  with  the 
current.  The  Auger  coefficient  is  extracted  as  a fitting  parameter 
in  the  numerical  solutions  that  simulate  the  experimental 

sublinearity. 

The  section  that  follows  presents  the  experimental  method, 
results  and  the  extraction  of  parameters  based  on  numerical 
solutions.  In  section  Section  5.3,  ve  discuss  the  sensitivity  and 
accuracy  of  the  method  and  provide  experimental  evidence  for  certain 
assumptions  used  in  the  interpretation  of  the  experimental  results. 
In  the  conclusion  section,  the  main  points  of  the  paper  are 
summarized. 

5*2  Experimental  Results  and  Numerical  Interpretation 
The  device  under  study  has  the  n+/p-/p+  structure  cut  from  a 
high  efficient  solar  cell.  The  p"  layer  has  the  doping  concentration 
Na  of  1.2  x 1015  cm"3  and  the  width  V of  100  pm.  The  device  is 
mounted  on  T0-5  package  and  bonded  with  16  bonding  wires  for  each 
electrode.  The  bonded  device  was  attached  to  the  heat  radiator  for 
the  possible  heat-up  of  the  device  during  the  high  injection  period. 
The  device  temperature  was  monitored  through  the  thermocouple 
mounted  on  the  heat  radiator.  We  applied  the  forward  bias  to  the 
diode  in  the  pulse  mode  to  prevent  it  from  heating  up,  as  shown  in 
Fig.  5.1.  The  repetition  rates  of  the  pulse  span  in  the  0.3  Hz  - 100 
Hz  range  depending  on  the  current  levels.  The  pulse  with  the  10~5 
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Fig.  5.1  Schematic  illustration  of  the  experimental  setup  for 
measurement  of  the  electroluminescence  versus  current 
for  p-i-n  diode.  MOS  switch  is  IRF510.  RL  is  0.2  E or  52  S 
depending  on  the  current  levels. 
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duty  cycle  of  3 sec  period  was  used  especially  at  the  high  current 
density  (J  > 1.5  x 1(P  A/ cm2).  The  electroluminescence  signal  was 
detected  by  InGaAs  photodetector.  The  detected  signal  was  amplified 
by  the  current  sensitive  amplifier  (PAR  184)  for  large  signals  or 
amplified  again  by  the  lock-in  amplifier  (PAR  124A)  for  the  weak 
signal  from  low  currents  (see  the  dashed  sequence  in  Fig.  5.1). 

Figure  5.2  shows  how  the  experimental  electroluminescence 
signal,  EL  as  a function  of  the  terminal  current,  I.  The  plot  has 
three  distinct  regions.  In  region  1,  EL  is  a superlinear  function  of 
I,  while  EL  depends  linearly  on  I in  region  2.  In  region  3,  the 
electroluminescence  signal  is  a sublinear  function  of  I and,  as 
discussed  later,  it  is  this  region  that  is  the  most  sensitive  to  the 
ambipolar  Auger  recombination  in  the  base.  To  understand  how  EL 
depends  on  I,  we  can  introduce  an  approximate  model  about  the 
carrier  transport  in  the  base  and  assume  that  the  carrier  plasma 
density  is  uniform  throughout  the  quasi-neutral  base.  Accurate 
numerical  results  that  allow  a more  precise  parameter  extraction 
will  be  presented  later  in  Section  5.2.2. 

5-2-l  Approximate  Analytical  Model 

In  this  approximate  model,  let  n be  the  uniform  electron-hole 
density  in  the  base  that  corresponds  to  the  terminal  current  I.  The 
density  of  electrons  and  holes  can  be  assumed  to  be  equal  for 
carrier  injection  levels  above  lO1*  cmA  Let  JQn  and  Jop  be  the 
saturation  current  densities  of  the  n+  and  p+  regions,  respectively, 
let  tsrh  be  the  Shockley-Read-Hall  (SRH)  ambipolar  base  lifetime  and 
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Fig.  5.2  Measured  electroluminescence  versus  current  through  the 
p-i-n  diode. 
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CA  be  the  ambipolar  Auger  coefficient  in  the  base.  For  carrier 
densities  high  enough  so  that  we  can  neglect  the  recombination  in 
the  space  charge  region  of  the  base  emitter  junction  as  well  as 
recombination  at  the  lateral  edge,  the  current  as  a function  of  n 
takes  the  form 

^ = A(don+JQp)n2/nj2  + eAWn/xgj^jj  + eAWCAn3,  (5.1) 

where  A is  the  area  of  the  device,  V is  width  of  the  p-  base,  and 

Jon  and  Jop  are  recombination  currents  in  heavily  doped  n+  and  p+, 
respectively. 

To  be  able  to  judge  the  relative  weight  of  the  three  terms  in 

(5.1),  we  can  use  values  for  the  recombination  parameters  that  are 

usually  adopted  in  the  literature.  The  sum  Jon  + JQp  is  of  the  order 
of  10-12  A/cm2 f the  high  injection  SRH  llfetime  is  usually  more  than 

100  ys  in  solar  cell  grade  material,  and  the  reported  values  of  CA 
are  in  the  range  from  0.4  x 10“30  to  2 x 10~30  cmV1  [49-54]. 

Therefore,  in  (5.1)  the  fist  term  dominates  when  n varies  between 
1016  and  1017  cm"3.  Thus,  in  this  range 

1 = A(Jon  + Jop)n2/ni2*  (5.2) 

At  the  same  time,  the  electroluminescence  signal  is  proportional  to 

the  integrated  pn  product  in  the  base.  If  F is  a proportionality 
constant  then 


EL  = Fn2. 


(5.3) 
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From  (5.2)  and  (5.3) 


EL  = F 


ni2I 


on  + 3op)A 


(5.4) 


and,  therefore,  in  the  range  from  1016  to  1017  cm"3  EL  is  a linear 

function  of  I.  This  range  roughly  corresponds  to  region  2 in  Fig. 
5.2. 

For  densities  above  1017  cm“3,  the  Auger  term  in  (5.1)  becomes 
important  and  now 


1 = A( Jon  + Jop)n2/ni2  + eAWCAn3 
1 = A(Jon  + J0p)EL/(Fni2)  + eAVCA(EL/F)3/2. 


(5.5  a) 
(5.5  b) 


Equation  (5.5  b)  implies  that  for  carrier  densities  above  1017  cm-3 
the  electroluminescence  increases  sublinearly  with  I and  this 
correspond  to  region  3 in  Fig.  5.2.  Region  1 corresponds  to  carrier 
densities  below  1016  cm"3,  where  (5.1)  becomes 


I-A(J0n  + 3op)n2/ni2  + eAVn/tSRH 

I=A(J0n  + Jop)EL/<Fni2)  + eAW(EL/F)1/2/xSRH 


(5.6  a) 


(5.6  b) 


which  implies  a superlinear  dependence  of  EL  on  I.  To  this 
superlinear  dependence  also  contributes  the  recombination  current  in 
the  space  charge  region  of  the  emitter  junction  as  well  as  the 
recombination  current  at  the  lateral  surface  of  the  diode.  This  is 
the  case  since  these  currents  depend  on  n according  to  na,  where  a 
is  in  the  vicinity  of  1. 
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The  above  arguments  would  be  more  convincing  if  we  could 
independently  calculate  the  average  carrier  density  as  a function  of 
the  current  and  verify  that  the  demarcation  levels  between  regions 
1,  2,  and  3 correspond  to  carrier  densities  in  the  vicinity  of  1016 
and  1017  cnr3.  Here,  by  average  carrier  density  we  mean  the  rms 
value  of  the  pn  product  in  the  base.  This  carrier  density,  n,  can  be 
calculated  from  (5.3)  provided  that  EL  can  be  measured  at  a bias 
level  where  n is  known.  Low  injection  can  serve  as  such  a bias  level 
since,  under  the  uniform  carrier  density  approximation, 

EL0  = Fno2  = Fni2exp(eV0/kT),  (5.7) 

where  VQ  is  the  terminal  voltage.  In  low  injection  the  current  is 
low  enough  to  eliminate  any  series  resistance  voltage  drop,  so  that 
(5.7)  holds.  In  Fig.  5.3,  we  plot  ELq  as  a function  of  the  terminal 
voltage  in  low  injection.  It  can  be  seen  that  EL0  changes  according 
to  (5.7)  and,  thus,  series  resistances  effects  are  negligible.  By 
combining  (5.3)  and  (5.7)  we  obtain 

n2  = (EL/EL0)ni2exp(eVo/kT).  (5>8) 

Therefore,  by  measuring  EL  at  an  injection  level  low  enough  so 
that  (5.7)  applies,  we  can  calculate  F.  Knowledge  of  F allows 
calibrating  the  measurements  at  high  current  levels  where,  due  to 
series  resistance  and  high  injection  effects,  the  terminal  voltage 
is  not  a useful  parameter.  From  Fig.  5.2  and  5.3  and  (5.8)  the 
boundary  between  regions  1 and  2 in  Fig.  5.2,  corresponds  to  n = 1.6 
x 1016  cm-3  while  the  boundary  between  regions  2 and  3 corresponds 


106 


TERMINAL  VOLTAGE  (mV) 


Fig.  5.3  Measured  electroluminescence  versus  terminal  voltage  of  the 
p-i-n  diode  at  low  current. 
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to  n = 9.8  x 1016  cm"3  in  agreement  with  the  previous  calculations 
based  on  generally  accepted  recombination  parameters.  Equation  (5.8) 
can,  in  addition  to  n,  provide  an  approximate  calculation  for  the 
sum  Jon  + JQp.  From  (5.2)  and  (5.8) 

3on  + 3op  = [ A(EL/ELg)exp(eV0/kT)  ] (5.9) 

where  EL  and  I belong  to  region  2.  From  (5.9)  and  Fig.  5.2  and  Fig. 

5.3,  Jon  + JQp  = 8.7  x 10-13  A/cm2,  in  agreement  with  the  original 
assumption  about  the  sum  JQn  + jQp. 

From  the  above  discussion,  and  since  Jon  + JQp  is  known,  one 
might  be  attempted  to  calculate  CA  from  (5.5)  and  (5.8).  However,  as 
we  point  out  later,  for  injection  levels  above  1017  cm"3  the  carrier 
density  is  not  uniform  and  the  above  approximations  can  lead  to 
large  errors.  In  order  to  derive  more  accurate  values  for  Jon,  JQp, 
and  CA,  we  will  simulate  the  experimental  results  using  numerical 
solutions  for  the  carrier  density  across  the  quasi  neutral  base.  The 
simulation  will  be  restricted  to  the  regions  2 and  3 because  these 
regions,  as  explained  above,  are  the  most  sensitive  to  Jon,  Jop  and 

CA,  and  they  are  much  more  immune  to  parasitic  currents  than  the 
region  1. 

5.2.2  Numerical  Analysis 

To  accurately  simulate  the  experimental  results,  the  transport 
equations  were  numerically  and  self-consistently  solved  from  the 

n+/p  to  the  p~/pf  edge  of  the  quasi  neutral  base.  The  transport 
equations  are 
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d*n 


Jn(x>  = -QUnn(x) 

dx 

(5.10  a) 

•Jp(x)  = qPpP(x) 

d't’p 

dx 

(5.10  b) 

p(x)  = Na  + n(x) 

(5.10  c) 

p(x)n(x)  = ni2exp[q(<f>p  - <f>n)/kT] 

(5.10  d) 

n(V)p(V) 

Jn(W)  = — 

n-iz 

^op 

(5.10  e) 

n(0)p(0) 

Jp(°)  = - Jon.  (5.10  f) 


The  solution  of  these  equation  provides  n(x)  and  I.  Since  n is  now  a 
function  of  x equation  (5.3)  becomes 


EL 


F 

W 


V 

n2(x)dx. 

Jo 


(5.11  a) 


As  we  prove  in  appendix  C,  to  include  the  free  carrier  absorption  we 
must  multiply  the  expression  in  (5.11  a)  by  following  factor. 

IV 


MF  = 1 - a 


n(x)dx. 


0 

f 

1 ) 

[V 

In 

+ y -ln(y)  - 

exp(-l/x)dx 

V J 

0 

a 

n(x)dx 

0 

(5.11  b) 
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where  a is  the  absorption  cross  section  of  the  electron  hole  plasma 
and  y is  arbitrary  number  much  lager  than  10.  The  free-carrier 
absorption  factor  is  written  here  in  terms  of  n,  but  this  is  of  no 
consequence  since  it  becomes  important  only  if  the  carrier  density 
is  above  1017  cm~3  where  n(x)  and  p(x)  are  essentially  the  same. 
Thus  (5.11  a)  becomes 


EL  = (F/V) 


MF 


' P 

n2(x)dx. 

JJo 


(5.11  c) 


As  in  the  previous  subsection,  the  proportionality  constant  F is 
calculated  from  the  voltage  dependence  of  EL  in  low  injection,  which 
is  shown  in  Fig.  5.3.  Instead  of  the  simplified  expression  in  (5.7), 
we  now  take  into  account  the  minority  carrier  variation  across  the 
base.  Since  the  base  is  uniform,  closed-form  expressions  are 
available  for  the  minority  carrier  density  in  the  base,  and  now 
(5.7)  becomes 


EL0  = Fn02  = Fni2exp(eV0/kT)  . 

(L/V)[sinh(V/L)  + (SL/D)cosh(V/L)  - SL/D 

u /tT  ,T  I ” • (5.11  d) 

cosh(W/L)  + (SL/D)sinh(W/L) 

In  (5.11  d),  S is  the  recombination  velocity  at  the  p-/p+ 
mterface  and  L is  the  low  injection  diffusion  length.  Since  in  high 
quality  solar  cells  SL/D  is  a small  number  and  L can  be  several 
hundred  microns,  the  large  fraction  in  (5.11  d)  is  very  close  to  1. 
Thus  (5.7)  overestimates  EL0  by  5%  to  10%  only. 
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5.2.3  Extraction  of  J„n  and  J^p 

To  reduce  the  number  of  unknowns  and  calculate  L and  S required 
by  (5.11  d)  for  the  calculation  of  F,  an  independent  experiment  was 
done  to  calculate  the  low  injection  lifetime,  x,  and  S.  The  diode 
was  exited  by  a uniformly  absorbed  pulse  of  a nanosecond  NYAG  laser 
( = 1.064  pm)  and  the  short  circuit  current  response  was  analyzed 
using  the  FDTA  method  [55].  The  amplitude  and  phase  response  of  the 
short  circuit  current  are  shown  in  Fig.  5.4  and  Fig.  5.5, 
respectively.  The  low  injection  lifetime  from  this  experiment  was  x 
= 50±5  ys  and  the  recombination  velocity  was  S = 26.5±5  cm/s.  From 
(5.11  d)  and  Fig.  5.3,  F = 4.63  x 10-29  (arb.  units).  The 
uncertainty  in  x and  S has  negligible  effects  on  F since,  as  already 
discussed,  the  big  fraction  in  (5.11  d)  is  close  to  1 and 
insensitive  to  L = 420  ym  and  SL/D  = 0.032.  Knowledge  of  S allows 
calculation  of  Jop  through  the  relation 

^op  = eSnj2/N^  = (4.4  + 0.9)  x 10~^  A/cm2.  (5.12) 

Since  Jop  is  independently  known,  fitting  the  experimental  results 
in  region  2 with  the  numerical  solutions  will  yield  JQn.  As  we 
demonstrate  in  the  next  section,  any  decrease  in  value  of  the 
mobilities  in  region  2,  has  no  effect  on  the  numerical  solutions  for 
the  EL  dependence  on  I,  in  the  same  region.  In  addition,  the  free- 
cairier  absorption  is  not  important  for  carrier  densities  below  1017 
cm"3'  Therefore,  Jon  is  essentially  the  only  fitting  parameter  in 
region  2 and  its  value  almost  uniquely  depends  on  the  experimental 
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Fig.  5.4  Amplitude  response  from  frequency-domain  transient  analysis 
for  n+/p"/p+  diode.  The  parameters  for  the  best  fit  are 
x = 50  ps,  S = 26.5  cm/sec. 
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Fig-  5.5  Phase  response  from  the  frequency-domain  transient  analysis 
for  n+/p“/p+  diode.  The  parameters  for  the  best  fit  are 
T = 50  ys,  S = 26.5  cm/sec. 
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results  and  Jop.  When  Jop  = 4.4  x 10  13  A/cm2,  fitting  in  region  2 
yields  Jon  = 4.3  x 10"13  A/cm2.  The  variation  of  JQp  in  (5.12) 
induces  a similar  variation  in  Jon.  The  effects  of  these  variations 
on  the  calculation  of  CA  will  be  discussed  in  the  next  section. 
Fitting  in  the  vicinity  between  regions  1 and  2 yields  the  high 
injection  SRH  ambipolar  lifetime  xSRH  = 150±20  ps.  The  uncertainty 
in  this  lifetime  has  negligible  effects  in  the  upper  half  of  region 
2 and  almost  no  effect  in  region  3. 

5.2.4  Extraction  of 

CA  is  obtained  from  the  fitting  in  region  3.  Here,  the  value  of 

the  mobilities  and  the  value  of  a do  affect  the  numerical  results. 

The  variation  of  mobilities  with  carrier  density  was  taken  from  the 

book  by  Baliga  [56]  while  the  value  for  a was  taken  from  data  by 

Svantesson  [57].  Svantesson  measured  a at  a wavelength  of  1.064  pm 

and  found  o = 5.1  x 10  ^ cm2.  The  electroluminescence  spectrum  of 

silicon  is  sharply  peaked  at  1.1  pm  and,  therefore,  we  can 

extrapolate  the  previous  value  of  a following  the  cubic  law 

dependence  of  a on  the  wavelength  [58].  Doing  so  yields  a = 5.64  x 
10~18  cm2> 

Using  these  values  for  the  mobilities  and  a and  having  already 
calculated  Jon  and  Jop,  the  fitting  in  region  3 becomes  a one 

parameter  fit,  CA  being  the  fitting  parameter.  CA  was  changed  by 
steps  of  0.01  x 10-30  cm6s-l.  The  best  fit  was  obtained  fof  ^ = 

1.77  x 10-30  cm^s-1.  Figure  5.6  shows  the  experimental  and  numerical 
results  for  a current  variation  over  three  orders  of  magnitude  that 
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Fig.  5.6  Curve  fitting  based  on  numerical  calculation. 

The  parameters  for  the  best  fit  are  Jon  = 4.3  x 10~13 
A/cm^.  JDp  * 4.4  x 10-13  A/cm^,  and  «=  1.77  x 10“30 
cm^s-l.  Mobilities  are  from  ref.  (56]. 
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covers  regions  2 and  3.  The  overall  agreement  is  better  than  9Q%. 
Figure  5.7  emphasizes  the  high  current  region  and  demonstrates  the 
sensitivity  of  the  fitting  to  CA.  In  Table  5.1  ve  summarize  the 
values  of  the  extracted  parameters. 


Table  5.1  Summary  of  the  data  set  for  the  best  fit. 


o 

3 

^op 

tSRH 

cA 

(10  13  A/cm^) 

(10-13  A/cm^) 

(US) 

Iq-30  cm6s-l 

4.3 

4.4 

150 

1.77 

5.3  Discussion 

5*3.1  Sensitivity  of  Ca  to  other  parameters 

Figure  5.8  demonstrates  the  effect  that  changes  in  the 
mobilities  have  on  the  numerical  calculations.  We  see  that  for 
currents  in  the  region  from  0.1  to  5 A,  the  numerical  results  are 
almost  completely  insensitive  to  any  changes  in  the  mobilities. 
Since  this  region  is  the  same  as  region  2,  we  can  safely  conclude 
that  the  extraction  of  Jon  does  not  depend  on  the  exact  knowledge  of 
the  mobilities. 

Figure  5.9  (a)  and  Fig.  5.9  (b)  show  that  the  plasma  densities 
in  the  p-  region  for  the  cases  of  with  and  without  the  Auger 
recombination,  respectively.  We  can  see  that  the  plasma  densities  at 
low  currents  are  same  for  both  cases.  However,  they  are  different  at 
high  currents.  The  plasma  distributions  without  the  Auger 
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Fig.  5.7  Sensitivity  of  fitting  to  Auger  coefficient. 

The  best  fit  is  for  = 1.77  x 10~30  cm^s-!* 
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Fig.  5.8  Sensitivity  of  fitting  to  mobilities. 

No  doping  variation  of  mobilities  is  assumed  for  the  upper 
curve.  The  bottom  curve  assumes  mobilities  from  ref.  [56]. 
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Fig.  5.9  Plasma  densities  across  the  p~  base  with  or  without  Auger 
recombination  for  different  currents. 

(a)  = 1.77  x 


(b)  = 0 cm^s 
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Fig.  5.9  (continued) 
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recombination  at  high  currents  is  almost  flat  through  the  p~  region. 
But  it  is  not  uniform  for  the  case  with  the  Auger  recombination.  The 
plasma  densities  for  this  case  at  x = 0 and  x = y are  larger  than  in 
the  middle  region  of  the  p~  layer,  because  the  ambipolar  diffusion 
length  becomes  shorter  due  to  the  Auger  recombination. 

Ve  also  checked  the  sensitivity  of  CA  to  the  mobilities  and  a 
in  this  experiment.  The  CA  varies  from  1.63  x 10-30  cm6s-l  to  1>88  x 
10-30  cm6s-l  £or  „ , 5.64  x 10_i8  cll-2±  20  Xi  1(  also  changes  1-50  x 

10-30  t0  1.94  x 10-30  cm6s-l  £or  , Mn  p ± 20  Xf  uhere  ^ are 
well  below  their  values  at  low  injection. 

The  effects  of  Jop  and  Jon  on  CA  are  summarized  in  Table  5.2. 

The  first  column  in  Table  5.2  represents  the  data  set  for  the  best 
fit.  The  second  column  shows  that  even  under  the  conditions  of  20  X 
higher  Jop  and  20  % less  Jon  than  the  values  in  the  first  column  the 
CA  is  same  as  in  the  best  fit  case.  The  third  column  is  for  the 
opposite  case  of  the  second.  In  this  case  we  have  the  slightly 
different  value  for  CA. 


5‘3*2  Band  Gap  Narrowing  and  Transition  Probabilities 

The  analysis  presented  in  section  Section  5.2  is  based  on  two 
implicit  assumptions.  Zero  band-gap  narrowing  and  an 
electroluminescence  intensity  that  is  proportional  to  the  pn 
product. 

We  checked  qualitatively  the  assumption  of  zero  band-gap 
narrowing  for  the  plasma  densities  under  study  by  measuring  EL  with 
and  without  interference  filter  between  the  p+/p-/n+  diode  and  the 
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Table  5.2  Sensitivity  of  the  saturation  current  to  fitting. 


(A)  Jop  (A/ cm2)  x 10-13 

Jon  (A/cmf)  4.3  x i0-13 


(cm^s- 

“0.121 

0.394 

0.816 

2.02 

7.6 

24.4 

46.0 

57.6 


) 1.77  x 10"30 

29779 
99519 

2.059  x 105 
4.995  x 105 
1.721  x 106 
4.614  x 10^ 
7.456  x 106 
8.753  x 106 
1.071  x 107 


ELo 

5.3  x 10~13 

3.3  x 10~13 
1.75  x 10~30 

29783 

99645 

2.062  x 105 
5.011  x 105 
1.732  x 106 
4.633  x 106 
7.469  x 106 
8.762  x 106 
1.071  x 106 


ELo 

3.5  x 10-33 
5.7  x 10~13 
1.96  x 10-30 

29316 

99151 

2.066  x 103 
5.044  x 105 

1.739  x 106 
4.633  x 106 
7.454  x 106 

8.739  x 106 

1.067  x 107 


77.2 
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sensor.  The  interference  filter  has  1.1  pm  transmission  center 
wavelength  and  sharp  cut-off  characteristics  for  other  wavelength. 
As  shown  in  Table  5.3,  the  ratios  of  two  ELs  at  various  currents  are 
almost  same.  This  suggests  that  band-gap  narrowing  in  the  range  of 

the  plasma  densities  under  study  is  either  not  existing  or 
negligible. 

The  second  assumption  is  equivalent  to  assuming  that  the 
probability  for  radiative  transitions  is  independent  of  the 
injection  level.  This  assumption  has  been  experimentally  verified  by 
Gerlach  and  et  al.  for  carrier  plasma  densities  up  to  2 x 101?  cm-3 
[59].  To  experimentally  verify  the  same  assumption  for  higher  plasma 
densities,  we  excited  a lightly  doped  silicon  wafer  with  the  second 
harmonic  (X  = 0.532  pm)  of  a nanosecond  pulsed  NdYAG  laser.  The 
energy  of  the  pulse  was  5 mJ  and  it  was  attenuated  by  a factor  of 
1000  by  a neutral  density  filter.  The  diameter  of  the  beam  was  3 mm. 
For  a reflection  coefficient  (=0.3)  and  an  absorption  coefficient  of 
0.9  x 103  cm"1,  right  after  the  pulse  the  plasma  density  at  the 
surface  becomes  about  1.2  x 10^  cm-3  and  decays  exponentailly 
inside  the  bulk  according  the  absorption  depth.  The 
electroluminescence  signal  right  after  the  pulse  comes  mostly  from 
the  excited  surface  and  the  first  micron  of  depth  where  the  plasma 
density  changes  from  1.2  to  0.45  x 10™  cm"3.  The  intensity  of  the 
laser  was  lowered  by  an  additional  factor  of  2 by  placing  another 
neutral  density  filter  between  the  first  and  the  wafer.  Thus  the 
detected  electroluminescence  comes  from  a plasma  density  varying 
between  1.2  and  0.22  x 10™  cm"3.  The  transient  electroluminescence 
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Table 


3.3  Electroluminescence  results  measured  with  (ELi)  and 

without  (EL2)  1.1  nm  filter.  EL^  and  EL2  were  measured  at 
different  sacles  of  lock— in  amplifier. 


CURRENT  (A) 

.EL!  (mV) 

(Without  filter) 

EL2  (mV) 

(With  filter) 

1.14 

22.2 

24.4 

5.4 

91.2 

101.6 

10.2 

158.4 

178.0 

40.0 

500.0 

556.5 

70.0 


744.0 


823.0 


LUMINESCENCE  (ARB.  UNIT) 
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USER  INTENSITY  (ARB.  UNIT) 


Fig.  5.10  Measured  luminescence  versus  laser  intensity. 
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signal  vas  detected  by  a InGaAs  detector  placed  from  the  other  side 
of  the  wafer  right  across  the  excited  spot.  This  is  a convenient 
configuration  because  it  completely  filters  the  powerful  laser  pulse 
through  the  wafer.  The  signal  was  amplified  by  a current  amplifier 
and  displayed  on  a scope.  Simultaneously,  the  laser  pulse  intensity 
was  monitored  by  employing  a beam  splitter  and  a Si  detector  for  a 
more  accurate  assessment  of  the  relative  intensity  of  the  exciting 
signal.  Figure  5.10  shows  the  peak  of  the  electroluminescence 
intensity  as  a function  of  the  laser  intensity.  The  figure  clearly 
demonstrates  that  the  electroluminescence  changes  according  to  the 
square  of  the  light  intensity.  Since  the  electroluminescence  signal 
comes  from  plasma  densities  varying  from  1.2  to  0.22  x 101®  Cm-3, 
and  the  carrier  density  immediately  following  the  pulse  is 
proportional  to  the  laser  intensity,  we  conclude  that  EL  is 

proportional  to  the  pn  product  for  carrier  densities  up  to  103-® 
cm~3. 


5.4  Conclusions 

The  recombination  lifetime  under  high  injection  conditions  in 
the  base  of  n+/p~/p+  diode  was  investigated  on  the  basis  of  the 
relation  between  electroluminescence  and  terminal  current.  For 
carrier  plasma  densities  up  to  1017  cm'3  the  relation  between  the 
electroluminescence  signal  and  the  current  through  the  device  is 
linear  and  dominated  by  the  heavily  doped  region  recombination.  This 
relation  becomes  sublinear  for  higher  carrier  density  due  to  Auger 
recombination  in  the  bulk  of  the  p-  base.  The  Auger  recombination 
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coefficient  is  extracted  by  fitting  the  experimental  sublinearity 
with  numerical  solutions  that  express  the  electroluminescence  signal 
as  a function  of  the  terminal  current.  In  the  1017  to  6 x 1017  cm"3 

carrier  plasma  density,  the  Auger  coefficient  was  found  to  be  1.77  x 
10~3®  cm^s-!. 


CHAPTER  SIX 

SUMMARY  AND  SUGGESTIONS 

This  dissertation  deals  with  the  device  and  physical  parameters 
of  bipolar  transistors.  The  device  parameters  are  dependent  on  the 
device  model.  The  parameters  investigated  in  this  work  are  for  the 
Gummel-Poon  model. 

The  base  resistance  was  obtained  by  measuring  the  input 
impedances  of  the  common-base  configuration  at  low  current.  The 
technique  includes  a correction  for  zeroing  errors.  The  correction 
is  based  on  comparing  the  directly  measured  common-base  input 
capacitance  with  that  obtained  by  combining  the  input  capacitances 
of  the  two  configurations  for  forward  and  inverse  operation.  The 
values  for  capacitances  measured  by  the  new  method  agree  with  the 
values  from  the  standard  method,  which  uses  accurate  dummy  devices 
for  zeroing.  The  advantage  of  the  new  method  is  that  the  accurate 
dummy  devices  are  not  necessary.  Furthermore,  the  accuracy  of  the 
new  method  is  about  +0.3  fF,  compared  to  1-2  fF  for  the  standard 
method.  The  standard  method  will  not  be  good  enough  for  future 
submicrometer-emitter  transistors  with  junction  capacitances  of  1-5 
fF.  The  junction  capacitances,  the  base  and  collector  saturation 
currents,  and  the  device  temperature  are  also  extracted  from  the 
detail  analysis  of  the  base  resistance  dependence  on  bias.  The 
method  is  also  demonstrated  on  a heterojunction  bipolar  transistor. 

A model  to  investigate  the  dependence  of  the  emitter 
lesistance  on  bias  and  the  structural  parameters  was  developed.  The 
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model  combines  the  emitter  and  contact  crowding  model  with  the 
distributed  model  for  a transistor.  The  emitter  resistance  was 
evaluated  by  the  circuit  simulator  SLICE  using  a power  dissipation 
method.  The  general  trends  of  the  emitter  resistance  agree  with  the 
experimental  results.  The  main  conclusion  of  this  topic  is  that  when 
the  emitter  crowding  becomes  significant,  the  emitter  resistance 
becomes  current  dependent.  For  the  simplification  of  the  analysis, 
the  current  gain  3 and  the  base  sheet  resistance  were  considered  to 
be  current  independent.  Their  dependences  can  be  included  in  the 
model  for  more  accurate  analysis. 

A method  for  extraction  of  the  collector  resistance  has  been 
developed.  The  basic  concept  of  the  method  is  to  use  the  substrate 
current  of  the  integrated  n-p-n  transistor.  This  current  is  either 
the  emitter  or  collector  current  of  the  parasitic  p-n-p  transistor, 
depending  on  biasing  conditions.  The  parasitic  transistor  is  turned 
on  by  the  voltage  drop  on  the  collector  resistance  of  the  n-p-n 
transistor.  The  separation  of  the  collector  resistance  components  is 
possible  by  operating  the  p-n-p  transistor  in  both  upward  and 
downward  modes.  It  is  difficult  to  accurately  estimate  the  value  for 
the  collector  resistance  at  high  current  because  it  is  strongly 
effected  by  the  high  injection  phenomena.  Since  the  new  method 
yields  the  current  dependence  of  the  collector  resistance,  it  may  be 
used  as  a tool  for  studying  the  high  injection  phenomena. 

The  recombination  lifetime  under  high  injection  condition  in 
the  base  of  n+/p~/p+  diode  was  investigated  on  the  basis  of  the 
relation  between  the  electroluminescence  and  terminal  current.  The 
Auger  coefficient  is  extracted  by  fitting  the  experimental  data  with 
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numerical  solutions  that  express  the  luminescence  signal  as  a 
function  of  the  terminal  current.  In  the  10^  to  6 x 10^^  cm-3 
carrier  plasma  density,  the  Auger  coefficient  was  found  to  be  1.77  x 
10  30  cm^s  1.  The  procedure  used  can  be  also  applied  to  a 
heterojunction  p/i/n  structure. 


APPENDIX  A 

INPUT  IMPEDANCE  FOR  COMMON-EMITTER  CONFIGURATION 
To  derive  Z^ne  we  use  the  simplified  equivalent  circuit  for  the 
CE  configuration  shown  in  Fig.  A. 2.  This  circuit  was  derived  from 
the  hybrid-rt  circuit  of  Fig.  A.l  by  neglecting  ry,  rc,  and  rQ.  We 
further  assume  that  rn»l /wCn.  The  validity  of  all  assumptions 
will  be  justified  below.  The  result  for  Zine,  neglecting  the 
higher  order  frequency  terms,  is 


'me 


f,  (1+&mre)cscK1+gmre)(2cy+CSC)+2cn]) 

rbx  + rbip  ~ — — 1 

l fclI+cy+Csc  +gmre(Csc+Cy)  ]2  j 

2 


+ r£ 


- (j/») 


(Cn+Cy+CSc)+gmre(cy+Csc)  > 

1+Smre 


(CJl+cy+csc)  + (Cy+csc)gmre. 


Assuming  further  that  gmre«l,  we  obtain 


(A.l) 


Jine  ~ rbx  + rbi 


cn+cy 


lCn+Cy+C$c. 


+rc 


Cji+Cy+CSci 


+ (l/ju) 


cn+cy+cSC 


(A. 2) 
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Fig. 


Fig. 


Csc 


A.l  Small-signal  equivalent  circuit  for  the  analysis. 


B rbx  VA  rbi  vB 


.2  Simplified  common-emitter  equivalent  circuit  obtained  from 
Fig.  A.l  by  neglecting  the  r^  and  r0,  and  the  collector 
series  resistance  rc. 


APPENDIX  B 

INPUT  IMPEDANCE  FOR  COMMON-BASE  CONFIGURATION  Zinb 
Using  the  simplified  diagram  of  Fig.  B.l  and  assuming 
rri>>l/wCTX,  rbx<<l/ioCgc,  ar*d  rbi<'<'l/wCp»  we  obtain 


zinb  “ 


(1+gmre)+jwCn(re+rbi+rbx) 
j “Cji+gm 


= re  + rbi 


( w2Crt2  ^ 


^m2+w2CII2 


Zr  2 


co^C 


+ r 


gm2+«2Cn2 


Sm 

gm2+to2Cn2 


^ &m(rbi+rbx) 

_j . (B.l) 

/wCn) 


Assuming  further  that  gra«wCrt,  and  gm(rbi+rbx)«l,  we  have 


Zinb  ~ re  + rbi  + rbx  ~ j(l/«Cn). 


(B.  2) 


The  validity  of  the  assumptions  used  to  derive  (A. 2)  and  (B.2) 
was  checked  by  calculating  the  input  impedances  Zine  and  Zinb  from 
the  complete  equivalent  circuit  shown  in  Fig.  A.l  using  SPICE.  The 
inputs  for  the  simulation  were  the  measured  parameters  rjj 
(determined  from  the  DC  IB-VBE  curve),  rbx,  rbi,  re,  Cn  and  Cjc. 
Reasonable  estimated  values  were  used  for  the  other  parameters, 
including  parasitic  capacitances.  The  agreement  between  the  simple 
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results  given  by  (A. 2)  and  (B.2)  and  the  circuit  simulation  vas 
better  than  1 percent  for  Vg£<0.45V  in  the  frequency  range  from  10 
to  100  MHz.  This  agreement  shows  the  excellent  accuracy  of  (2.1)- 
(2.4)  used  in  this  work. 


Fig.  B.l  Simplified  common-base  equivalent  circuit. 


APPENDIX  C 

FREE  CARRIER  ABSORPTION 


The  average  number  of  photons  emitted  under  an  angle  9 that  is 
not  absorbed  by  free  carriers  is  proportional  to 


exp(-<r 


V 

n(y) 

0 


dy 

d9) 

cos9 


(C.l) 


where  a is  absorption  cross-section,  n(y)  is  carrier  density,  and 
others  are  defined  in  Fig.  C.l. 

The  total  number  of  photons  that  are  emitted  under  this  angle 
and  leaves  the  cell  surviving  the  absorption  is  proportional  to 


exp(-crn/cos9).2nsin9d9 


where  crn 


a 


V 

n(y)  dy. 

0 


(C.2) 

(C.3) 


The  2nsin9d9  in  (C.2)  is  a geometrical  factor  in  the  spherical 
coordinate  for  taking  into  account  the  light  injected  into  all 
directions.  The  assumption  used  in  the  above  derivation  is  the 
total  attenuation  is  much  smaller  than  1. 

Then,  the  total  flux  leaving  the  cell  is  proportional  to 


'n/2 

exp(-CTn/cos9)2nsin9d9. 

JO 


(C.  4) 
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FRONTSIDE  (AR  COATED  LAYER) 


BACKSIDE  (SILVER) 


Fig.  C.l  Schematic  diagram  for  the  analysis  of  free  carrier 
absorption . 
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After  changing  variable,  (C.4)  becomes 


p/"n 

exp(-l/x)dx. 

0 


(C.  5) 


For  a certain  number  A which  is  larger  than  1,  (C.4)  is  expressed 
as 


t 

'A 

1/^n  ] 

an 

exp(-l/x)dx  + 

exp(-l/x)dx 

< V 

0 

A J 

'A 

P/ffn 

exp(-l/x)dx  + 

(1  - l/x)dx 

0 

A J 

'A 

In  (l/an)  + A -In  A - 

exp(-l/x)dx 

0 J 

(C.  6) 


The  measured  electroluminescence  is 


/ 

t 

'A  y 

1 - CTn 

In  (l/crn)  + A -In  A - 

exp(-l/x)dx 

\ 

0 JJ 

(C.  7 ) 


where  EL  is  calculated  electroluminescence  without  the  free  carrier 
absorption. 
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